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Preface 



Quantum process tomography is one of the most fundamental tool for 
the experimental study of open dynamics of a quantum system. Moreover, 
quantum process tomography is an important practical tool for implement- 
ing quantum information processing devices. By the very nature of quantum 
mechanics, quantum information processing cannot be achieved in isolation. 
The processing device must interact with a measuring apparatus to yield in- 
formation. Furthermore, the processing device, in general, is also affected by 
its surroundings (an environment) in an uncontrollable manner. This results 
in the necessity of treating the dynamics of such devices with the quantum 
theory of open systems. 

The basic building blocks of quantum processing devices have only been 
experimentally realized in the last ten years; hence quantum process tomog- 
raphy is a relatively new procedure. On the contrary, the theoretical analog 
of quantum process tomography is the dynamical map formalism, introduced 
almost fifty years ago by Sudarshan, Mathews, and Rau. During the last five 
decades this formalism has matured significantly. In some of the recent studies 
of the dynamical maps, a significant emphasis is placed on the reduced dy- 
namics of initially correlated states. While for quantum process tomography 
it is usually assumed that at the beginning of the experiment the state of the 
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system is uncorrelated with its surroundings. Based on some of the theoretical 
results, we study quantum process tomography for initially correlated systems. 

Since quantum process tomography is an experimental procedure, we 
are forced to deal with the idea of the preparation of initial states. In fact, one 
must analyze the role of the preparation procedures for any quantum experi- 
ment that interacts with its surroundings non-trivially. There are two different 
ways by which the preparation procedure can play a significant role in an open 
quantum experiment. First is due to the initial correlations between the states 
of the system and the environment. While the second deals the consistency of 
the preparation procedures; even when there are no initial correlations between 
the system and the environment. 

The starting point of this dissertation is a review of the dynamical 
map formalism and several quantum process tomography procedures. At this 
point, the mathematical structure of preparation procedures is discussed in 
detail. This allows us to investigate the role of preparation in quantum process 
tomography without any assumptions. 

The study of quantum process tomography with preparation proce- 
dure leads to another process tomography method that is independent of the 
preparation procedure. Furthermore, this procedure leads to a surprising re- 
sult; the map arising from this procedure (we call M-map) leads a quantitative 
measure for the non-Markovian memory effect on the system due to the ini- 
tial correlations between the system and the environment. A measure for the 
non-Markovian memory effect can play a crucial role in a coherence control 
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scheme. 

Many of the results deahng with negativity and non-hnearity in quan- 
tum process tomography are presented by concrete examples, rather than by 
a general mathematical analysis. Based on the results of these examples we 
try to generalize the role of preparation in quantum process tomography. On 
the other hand, the theory of preparation procedure, development of M-map 
and the memory matrix are worked out by general mathematical analysis. 

Kavan Modi, Austin Texas. 
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In recent years there has been a significant development of the dynam- 
ical map formalism for initially correlated states of a system and its environ- 
ment. Based on some of these results, we study quantum process tomography 
for initially correlated states of the system and the environment. This is be- 
yond the usual assumption that the state of the system and the environment 
are initially uncorrelated. Since quantum process tomography is an experi- 
mental procedure, we wind up having to study the role of preparation of input 
states for open quantum experiments. We work out a theory for the gen- 
eral preparation procedure, and study two preparation procedures in detail. 
In specific, we study the stochastic preparation procedure and the projective 
preparation procedure and apply them to quantum process tomography. The 
two preparation procedures describe the ways to uncorrelate the state of the 
system and the environment. However the specifics of how this is implemented 
plays a role on the outcomes of the experiment. 

X 



When the stochastic preparation procedure is apphed properly, quan- 
tum process tomography yields a linear process maps. We point out what it 
means to apply the stochastic preparation procedure properly by constructing 
several simple examples where inconsistencies in preparations leads to errors. 
When the projective preparation procedure is applied, quantum process to- 
mography leads to a non-linear process map. We show that these processes 
can only be consistently described by a general dynamical map, which we call 
M-map. The M-map contains all of the dynamical information for the state 
of the system without the affects of a preparation procedure. By carefully 
extracting some of this dynamical information, we construct a quantitative 
measure for the memory effect due to the initial correlations with the environ- 
ment. 



xi 



Table of Contents 



Acknowledgments |v] 

Preface Iviil 

Abstract E] 

List of Figures Ixvl 

Chapter 1. Introduction [1] 

1.1 Organization of this dissertation |3] 

Chapter 2. Quantum theory of open system [6] 

2.1 Closed dynamics IE] 

2.2 Open dynamics [TU] 

2.3 Dynamical map formalism [13] 

2.3.1 yi-form M 

2.3.2 S-form M 

2.3.3 Semi-group property 

2.3.4 Positivity classes [18] 

2.3.5 Choi representation [19] 

2.4 Dynamical maps from contraction [20] 

2.4.1 Initially product states [21] 

2.4.2 Initially correlated states [22] 

2.5 Size of the environment [29] 

2.6 Discussion [30] 



xii 



Chapter 3. Quantum process tomography 1321 

3.1 Quantum process tomography: general description [33] 

3.2 Standard quantum process tomography [35] 

3.3 AnciUa assisted quantum process tomography [37] 

3.4 Direct characterization and other methods [IQ] 

3.5 Discussion [H] 

Chapter 4. Preparation of input states 1431 

4.1 Preparations in quantum mechanics [S] 

4.1.1 Trace preservation [15] 

4.2 Preparations in open quantum mechanics [16] 

4.3 Two common preparation procedures [IS] 

4.3.1 Stochastic preparation [H] 

4.3.1.1 For open systems [IH] 

4.3.2 Projective preparation [51] 

4.3.2.1 For open systems [52] 

4.4 Discussion [55] 

Chapter 5. Preparation and quantum process tomography 1561 

5.1 Quantum process tomography experiment in steps [57] 

5.2 Stochastic preparation [5H] 

5.2.1 An example with stochastic preparation [32] 

5.2.2 An example with multiple stochastic preparations .... [60] 

5.3 Projective preparation [M] 

5.3.1 An example with projective preparation [66] 

5.4 Discussion [HH] 

Chapter 6. More on stochastic preparations [71] 

6.1 Pseudo-pure states [75] 

6.2 Negative maps due to control errors [75] 

6.3 Discussion [77] 



xiii 



Chapter 7. Dynamical M-map 1791 

7.1 Uninterrupted dynamics [H2] 

7.2 Memory due to correlations [83] 

7.3 An example with M-map and the memory matrix [86] 

7.4 Discussion [88] 

Chapter 8. Analysis of experiments 1891 

8.1 Experiment by Myrskog et al [89] 

8.2 Experiment by Howard et al [9T] 

Chapter 9. Conclusion and future directions [94] 

9.1 Future directions [M] 

Appendix 1961 

Appendix A. Quantum state tomography 1971 

Appendix B. M-map process tomography 1981 

B.l The qubit case M 

B.2 Beyond one qubit [T02] 

Appendix C. Experimental recipe for M-map 11031 

Appendix D. Notation [TU51 

Bibliography IllOl 

Index [m 

Vita 

xiv 



List of Figures 



2.1 The total state of the system and environment (p^^), with the 
state of the system (p^) is represented by color green and the 
state of the environment (p^) is represented by color red. The 
fuzzy part in-between represent the correlations between the 
system and the environment (p^^ — p^ (8) p^). The dotted blue 
lines represent the set all physical states in that space. Due to 
the correlations not all possible states of the system and the 
environment are allowed. Only the states that are compatible 
with the correlations are allowed, represented by the green and 

red ellipses on the bottom [TT] 

2.2 [Top] The total state p^^ evolves unitarily, which changes the 
polarization of the system, the environment, and correlations 
between them (note the color changes). [Bottom] While the 
reduced state of the system does not evolve unitarily. Note 
that the shape and the size of the set of states of the system 
(green ellipses) change from the initial time to the final time. 

We want to be able to describe this transformation [13] 



2.3 The eigenvalues of the dynamical map in Eq 2.43 are plotted 
as function of 2ujt. One of the eigenvalue is negative for certain 
values of a;t; we have taken C23 = 0.5. The negative eigenvalue 
is due to the initial correlations between the system and the 
environment [22] 

2.4 [Left] The total state of the system and the environment has 
correlations, therefore the set states that are physically available 
for the system part (green ellipse on bottom) are less than set of 
all state (dotted blue circle). This is the compatibility domain. 
[Right] Conversely, if one demands to embed set of all system 
states into the total state (compatible are in the green ellipse 
and incompatible are represented by the blue part of the circle 
on the bottom), then the total state may become unphysical 
(fuzzy blue and red area above) [27] 
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2.5 [Top] The total state evolves unitarily, which changes the po- 
larization of the system, the environment, and the correlations 
between them (note the color difference). [Bottom] The dynam- 
ical map acting on the compatibility domain results in physical 
states (green ellipses from left to right). But its action on set 
of all states (dotted blue circle on left) can take physical states 
to non-physical states (dotted red ellipse on right) EH] 

4.1 At the beginning of the experiment an unknown state of the 
system and the environment is present. The mth state of the 
system is then prepared, which diminishes the correlations be- 
tween the system and the environment HT] 

4.2 The figure above shows how that for two different projective 
preparations lead to the desired input states, but the corre- 
sponding states of the environment can be different. In both 
cases the initial total state is the same [Ml 



5.1 The eigenvalues of the process map in Eq. |5.8| are plotted as 
function of 2ujt. As expected the eigenvalues are always positive. 
The stochastic preparation procedure allows the experimentalist 
to prepare any pure state for the system. By convexity, then all 
possible states of the system are in the compatibility domain of 



the process map in Eq. 5.8 EH 



5.2 The eigenvalues of the process map in Eq. 5. 15| plotted as func 



tion of 2ujt. One of the eigenvalue is negative for certain values 
of ut. The negativity is due to the inconsistency in the prepa- 
ration procedure, and not due to the initial correlations with 
the environment 



5.3 The eigenvalues of the process map in Eq. |5.26| are plotted as 
function of 2ujt. One of the eigenvalue is negative for certain 
values of ujt\ we have taken = 0.5. The negativity here is 
due to the initial correlations and the inconsistence preparation 
procedure 

6.1 The eigenvalues of the process map found using the output 



states in Eq. |6.15| and the duals in Eq. |3.8[ One of the eigen- 
value is negative for certain times; we have taken e = 0.1. The 
negativity is due to the errors in the unitary operation imple- 
mented to prepare one of the input states [77] 
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Chapter 1 
Introduction 



Quantum information processing promises powerful computational meth- 
ods that surpass the methods of classical information processing PQ |2] . These 
methods rely on taking advantage of quantum parallelism by using quantum 
superposition and quantum entanglement as resources. In order to implement 
such a device one must have precise control over the system, and isolate it 
from the surrounding environment to preserve coherence. For delicate systems 
of this sort, it is nearly impossible to isolate the system of interest completely 
from its surroundings, while having a great deal of control. 

With the rising interest in quantum computation and quantum informa- 
tion processing, quantum coherence experiments are performed readily these 
days, though with relatively small systems. One of the major problem with 
these experiment is the loss of coherence due to the interaction between the 
system of interest and the unknown environmental states. The methods for 
studying the interaction between the system and the environment are given 
by the quantum theory of open systems. 

The quantum theory of open systems got its start in almost fifty years 
ago with the introduction of dynamical maps [U H] due to Sudarshan, Math- 
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ews, and Rau. Fifty years after its conception, the dynamical map formalism 
is finally being tested in the laboratory setting. The experimental determina- 
tion of a dynamical map is achieved by a procedure called quantum process 
tomography 

Any experiment, including quantum process tomography experiments, 
requires a method to prepare the initial states of the system at the beginning 
of the experiment [5]. We study the affects of the preparation procedure on 
quantum systems that interact with an environment The act of preparation 
has been neglected from the theory of quantum process tomography (and for 
all quantum experiments that interact with a non-trivial environment). We 
investigate this issue for quantum process tomography in detail in this dis- 
sertation. We present several simple examples to motivate our arguments. 
Based these arguments, we analyze two recent quantum process tomography 
experiments and show their results to be consistent with our analysis. 

In the course of our studies of the role of preparation procedures in 
quantum process tomography, we derive a new powerful method of quantum 
process tomography that is not affected by the preparation procedure. An 
experimental recipe to carry out this procedure is given at the end. 

The added advantage to this new procedure is an expression quantifying 

^Tomography, in the traditional sense, means to determine the internal structure of an 
object. Here the process is thought to be a physical object that interacts with the system 
for an instant, hence process tomography means the action of the object onto the system. 

^An environment is any degree of freedom that develops in time with the state of the 
system. 
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the memory due the initial correlations in the dynamics of the state of the 
system. Quantifying the memory is the crucial step to develop a coherence 
control scheme. If one can separate the purely dissipative terms from the terms 
that periodically recur, then a scheme can be developed to make use of the 
coherence periodicity, in battling decoherence. 

1.1 Organization of this dissertation 

We start with a brief review of closed dynamics of quantum systems 
in Chap. |2} We then motivate the necessity for studying open quantum 
dynamics. In this dissertation we work in the dynamical map formalism to 
study the open dynamics of quantum system^ Next we define the positivity 
classes associated with the dynamical maps and present an example of a not- 
completely positive dynamical map. The not-completely positive nature is 
attributed to the initial correlations between the system and the environment. 
We use this example later on as an inspiration in analyzing quantum process 
tomography experiments. 

In Chap. [3] we review several existing methods of quantum process 
tomography. We discuss the reasons for the different methods of quantum 
process tomography by analyzing the pluses and the minuses associated with 
each. We point out the central assumption in all of these procedure; that, 
the state of the system and the environment is uncorrelated at the beginning 

^See [HI [7] for other approaches to quantum theory of open systems. 
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of the experiment. When this assumption does not hold, a quantum process 
tomography experiment can yield a nonsensical process map for some cases. 

In Chap. |4] we present one of the central points of this dissertation. 
We discuss the general theory for preparing input states in terms of stochastic 
maps. We discuss in detail two of the most common preparation procedures 
practiced and compare them for closed and open quantum systems. For open 
systems, if care in not take in implementing a preparation procedure, then the 
state of the environment will pickup a non-trivial dependence on either the 
prepared state or the preparation procedure itself. This dependence can lead 
to nonsensical experimental results. 

In Chap. |5] we revisit the quantum process tomography procedures 
armed with the preparation techniques. We analyze the role that preparation 
procedure plays in a quantum process tomography experiment. We show that 
for certain types of preparation procedures, the quantum process tomography 
will fail to obtain a process map that correctly describes the physical process. 
We present three detailed examples dealing with these difficulties. 

In Chap. [6} we analyze the case when the preparation procedure itself 
is poor. We present several more examples with different scenarios of poor 
preparation procedures. The causes for the errors in this chapter are different 
from the causes discussed in Chap. |8j We compile a list of operations that 
can lead to errors for a quantum process tomography experiment. 

We offer a resolution by introducing a new process map that we call 
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dynamical M-map in Chap. [7j To obtain this map we extract the prepara- 
tion procedure from the dynamics. In this fashion the M-map contains all 
of the dynamical information for the system unaffected by the preparation 
procedure. We further show that, the M-map contains information about the 
dynamics of the system when it is correlated and when it is uncorrelated with 
the environment. This allows us isolate the reduced dynamics of the initial 
correlations between the system and the environment known as the memory 
due to the correlations and construct a quantitative measure for the mem- 
ory. Finally we give an example of M-map, and from it calculate some of the 
underlying dynamics. 

In Chap. [8} we analyze two experiments where the preparation proce- 
dures were not carried out consistently. The negative results obtained in these 
experiments are foreseen by our theory. We give the concluding remarks and 
potential future directions in Chap. |9| 

In App. |A| we briefly discuss quantum state tomography. In App. |B] 
we lay out a procedure to experimentally determine the M-map, and lay out 
a recipe for an experiment for a qubit system in AppjCj Finally in App. [D]we 
lay out all of the notation used throughout this dissertation. 
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Chapter 2 



Quantum theory of open system 



Pure states in quantum mechanics are generally represented by rays 
in the Hilbert space. Rays adequately explain the phenomena of quantum 
superposition, but they cannon represent classical mixtures of quantum states. 
For instance, consider a beam of particles with 70% in state |1) and 30% 
in state |0); this situation is not the same as quantum superposition. We 
can represent the state of the beam by a density matrix as p = 0.7 |1) (1| + 
0.3 |0) (0|. Here 0.7 and 0.3 represent the classical probabilities of finding the 
system in state |1) or |0). Any state of a quantum system can be written in 
terms of a density matrix^ be it a pure quantum state or a classically mixed 
one. 

Since density matrices describe classical mixtures of quantum states, 
they must have the following three properties. 



The eigenvalues of the density matrix represent classical probabilities for each 
pure state given by it's corresponding eigen-projectors. Thus, the trace con- 

'^We only consider finite dimensional systems throughout this dissertation. 



Tr[p] = 1 



Normalization, 

Hermitian, 

Positivity. 



P = P' 
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dition gives us the conservation of probabilities, Hermiticity of p guarantees 
that the eigenvalues are real, and finally with the last condition we demand 
that the eigenvalues are positive. 

Density matrices form a convex set; if {p(^)} are a set of density matrices 

then 

p = ^py^) (2.1) 
j 

is also a density matrix as long as pj are real, positive, and YljPj ~ 1- ^'^y 
density matrix that cannot be written in terms of a convex sum is called 
extremal. Extremal states are also pure states. An important decomposition 
is the eigen-decomposition, 

P = E^^-|^^-)<^^-|' (2.2) 

j 

where the Xj and {ipjl are the eigenvalues and eigenprojectors of p respec- 
tively. 

For a qubit (a two-level system), we can write the density matrix as 
the following 

p = ^{I + aidi + a2a2 + 030-3}, 
where I is the 2x2 identity matrix, aj are the Pauli spin matrices given by 

^^=(1 0)' ^^=(. ^=^=(0 -1)' 

and (01,02,03) are the components of the Bloch vector [8J satisfying \a\ < 1. 
The Bloch vector a has a nice geometric interpretation for a qubit; it can be 
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thought of as a vector inside of the unit sphere. The magnitude of a represents 
the polarization of the density matrix. For instance, if |a| = 1 then p is a pure 
state and if |a| =0 then p is a completely mixed state, with everything else 
in-between. 

In matrix form the density matrix takes the form 



with eigenvalues — 

When a state of the system evolves with an environment, it exchanges 
physical quantities such as polarization or phase with the environment. It 
also becomes correlated with the state of the environment. Such dynamics 
cannot be described adequately by a pure state; therefore density matrices are 
employed to study open dynamics. 

2.1 Closed dynamics 

The dynamics of a density matrix is governed by the von-Neumann 
equation 



where H{t) is the time dependent Hamiltoniarj^ By integrating the above 
equation, we find that the time evolution is also be expressed by a unitary 

^We take h = 1 from here on. 




1 / 1+0,3 CLi — ia2 

2 I fli + m2 1 - as 




ih 



dp 

di 



[Hit),p], 



(2.3) 
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operator, [/, as follows 



p{t) = U{p{to)) = U{t,to)p{to)U^it,to). 



(2.4) 



The unitary operator, U{t, to), takes the density matrix from time to to i in a 
linear fashion 



where ^jPj — 1- That is, the unitary operations preserve the mixture. An 
unitary operator in general can be written as 



where ^ represents time ordering and H{t') is the time dependent Hamilto- 
nian. In this dissertation we will only work with time independent Hamilto- 
nians for simplicity, though all results apply to time dependent cases as well. 

In quantum information theory, the information is processed by passing 
quantum states through various quantum logic gates, which are represented by 
unitary operations. Then we are interested only in the finite time evolution of 
a quantum state. Therefore, it is far more convenient to describe the evolution 
of quantum states with unitary operators rather than with the von Neumann 
differential equation. 

In the next section we investigate the characteristics of an open quan- 
tum system. Our goal is to describe the most general time evolution of the 
system of interest. The most general evolution includes non-dissipative evolu- 
tion, dissipative evolution, and the combination of the two. 




(2.5) 




(2.6) 
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2.2 Open dynamics 

In general, a state going through a quantum logic gate is not only sub- 
ject to the gate, but also to the surrounding environment. Such unwanted 
influences lead to noisy information and processing errors. This noise should 
be minimized for optimal performance. To optimize a gate we must first un- 
derstand the additional infiuences from the environment on the state. The 
presence of an interacting environment leads the system to experience "open 
dynamics", where physical quantities such as polarization and phase are af- 
fected by the state of the environment. 

More specifically, a quantum state in dynamics can experience changes 
in the relative phase between the orthogonal components as well as changes 
in its polarization. We can recognize that the magnitude of the polarization 
of the system density matrix alone does not change under unitary transfor- 
mations; that is eigenvalues of the density matrix are invariant under unitary 
transformations but not the eigenvectors. 

For a qubit, the Bloch vector under a unitary transformation goes 
through a rotation in the Bloch sphere, a(to) a'{t). While the magnitude of 
the Bloch vector remains unchanged, |a(to)| = In open evolution, the 

Bloch vector will experience rotations as before; additionally its magnitude 



may change in time, |a(to)| 7^ \(^'{t)\ (see Fig 2.2 for a graphical explanation). 

We can write down the interaction between the system and environment 
easily by treating the combined evolution in the closed from. Consider a 
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Figure 2.1: The total state of the system and environment (p^^), with the 
state of the system (p^) is represented by color green and the state of the envi- 
ronment (p^) is represented by color red. The fuzzy part in-between represent 
the correlations between the system and the environment (p^^ — p® ® p^). The 
dotted blue lines represent the set all physical states in that space. Due to 
the correlations not all possible states of the system and the environment are 
allowed. Only the states that are compatible with the correlations are allowed, 
represented by the green and red ellipses on the bottom . 



bipartite state p^^ of the system (labeled by §) and the environment (labeled 
by £). The total unitary evolution is as follows 

pS£(t) = Up''{to)Ul (2.7) 

For generality, the state of the system and the environment is assumed to be 



11 



correlated initially (see Fig. 2.1). We also assume that we cannot observe 
the state of the environment and do not know what the interaction unitary 
transformations are. If we had the knowledge of the state of the environment 
and the unitary transformations, by calculating the closed evolution we would 
know what the state of system will be at any point. Our only knowledge 
is of the initial and final states of the system. The initial and final states 
can be obtained by averaging over the environment in the above equation. 
Mathematically this corresponds to a partial trace operation 

P%) = Tt,[p'%)]. (2.8) 



See Fig 2.2 for a graphical explanation. 



There is no general equation in the differential form (like Eq. 2.3) that 
governs the dynamics of the state of the system from to to ^ It is easy to 
see that the mixture of the states of the system initially will not be the final 
mixture. For instance, p^(to) could begin as a pure state, and evolve into a 
mixed state, or vice versa. This is due the exchanges of physical quantities, 
such as polarization, between the system and the environment. They may also 
become correlated by loosing polarization to correlations. This exchange can 
be periodic^ semi-periodic, or purely relaxing. For periodic and semi-periodic 



■^If the process is Markovian in nature then the Kossakowski-Lindblad equation can be 
appHed O [TOj [HI [12] . A Markovian process is when the state of the system only depends 
on the state in the previous step. 

^Systems that have a memory are called non-Markovian. In this case the memory effect 
is due to the correlations with the environment [15] . 
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Figure 2.2: [Top] The total state p^^ evolves unitarily, which changes the 
polarization of the system, the environment, and correlations between them 
(note the color changes). [Bottom] While the reduced state of the system does 
not evolve unitarily. Note that the shape and the size of the set of states of 
the system (green ellipses) change from the initial time to the final time. We 
want to be able to describe this transformation. 



cases, the system may start pure, become mixed at some intermediary step, 
and then become pure or almost pure once again. 

2.3 Dynamical map formalism 

Suppose we are not interested in the dynamics of the system for all times 
but rather in the transformation of the state from time to to time t. Then we 
need to define a mapping from density matrices to density matrices such that 
all allowed initial states of the system are mapped to the corresponding final 
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states in a linear fashion. We closely follow the arguments originally put forth 
by Sudarshan et al. [3]. 



2.3.1 ^-form 

Consider an operator acting on the density matrix mapping the state 
to another density matrix linearly 



Above is labeled by two indices and the operator A is labeled by four, 
meaning ii is a. d x d matrix then ^1 is c?^ x c?^ matrix. We can think of 
the above equation as a super matrix A acting on a column vector p^. This is 
very similar to a classical stochastic process, where a stochastic matrix maps 
a probability vector to another probability vector. Matrix is A called the 
stochastic map [31 [H] . 

The only restriction we need to place on the stochastic map is that it 
map a density matrix to another density matrix. This implies that it must pre- 
serve trace, Hermiticity, and positivity of the density matrix. These restriction 
translate into the following properties for A 




(2.9) 




> 



Trace preservation, 
Hermiticity preservation, 
Positivity. 
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2.3.2 3-form 



Following Sudarshan et al. [3] again, let us define a dynamical map 
from the stochastic map by 



(2.10) 



For a 4 X 4 map, matrix A tranforms as follows 

/ All Ai2 Ays Ai4 \ ( All Ai2 A21 A22 \ 

Ais Ai4 A23 A24 



A21 A22 A23 A24 
A31 As2 Ass A34 
y A41 A42 A43 A44 j 



Azi Az2 A41 A42 
\ Azz Az4 A4Z A44 I 



(2.11) 



The properties of the stochastic map translate in the following manner 
for the dynamical map [3]. 



(J^ss' ,rr' ) 



Trace preservation, 
Hermiticity preservation, 
Positivity. 



The advantage of writing the dynamical map B is that it has nicer 
properties than the stochastic map A. The dynamical map is Hermitian and 
its diagonal d x d block elements have unit trace. For these nicer properties 
we have sacrificed the simplicity of the composition of the stochastic map on 
the state. The stochastic map acts as a matrix on a state that is a column 
vector, while the dynamical map acts on the state in the following manner 

'Brr';ss'pf,Ato)=pfsit)- (2-12) 

The consequence of sacrificing the simple action of the stochastic map on 
density matrices is that we have sacrificed the semi-group property that comes 
with it. We will discuss this matter in the next section. 
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First, we should remark that since the dynamical map can transform 
all possible states to all possible states, it describes the most general evolution 
of a quantum state for finite time. We need nothing more to to describe any 
time evolution of a state. Furthermore, the set of density matrices form a 
convex set, then the set of trace preserving positive maps that act on density 
matrices also form a convex set. Then any map written as 

^ = E^^-^^- ; o<^^- ; E^^^i (2.13) 

is a valid trace preserving positive map, given that all !Bj positive and preserve 
trace. Any positive maps that cannot be written as a convex sum of other 
positive map is called an extremal map. 

2.3.3 Semi-group property 

Suppose map A'^^^ti^to) takes a state from to to t and map A^'^\t2.,t) 
takes the state from t\ to ^2- Then the map that takes the state from to to t2 
is simply 

Plsit2) = A%,^,{hM)pl'Ati) (2-14) 

= 42U'(^2, tM"Ah) (2.15) 

= 4t"."(^2,to)p.V(to), (2.16) 

where A'^^^ = ^^^•'^1*^^^ simple by matrix multiplication. 

While for the corresponding dynamical maps we have 

pfM = Sf;.,,,(i2,ii)pfv(^i) (2.17) 
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= ^fXsAi2,h) o ^%,,,,Atuto)p%Ato) (2.18) 

= ^^XAt2,to)pf^>Ato), (2.19) 

where o represents composition of two maps, and not simple matrix multipli- 
cation. The same simple composition property of the stochastic maps does 
not hold for the dynamical maps, i.e. 

3(3) _^ 3(2)3(1) _ In general, knowing 
the mapping from to to t does not tell us very much about the dynamics from 
t to t2- This is due to the fact, whatever the correlations that the system and 
environment build up during their interaction from time to to t can comeback 
and play a role in the dynamics of the system from time t to t2- This memory 
effect in the dynamics results in non-Markovian dynamics [HI [13] . 

A special but important case is when the dynamical map has the semi- 
group property. The resulting dynamics is Markovian. A Markovian system 
has no memory, meaning the state of the environment and the correlations 
with the environment do not affect the dynamics of the system. In that case 
we can simply write down the state of the system at tn as 

p^itn) = B(t„,t„_i)o...oS(t2,t)oS(ti,to)p^(to) (2.20) 
= m„to)rp^to). (2.21) 

This case leads to the well known Kossakowski-Lindblad master equation [Qj 

nnj. 
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2.3.4 Positivity classes 

One advantage of the dynamical map is that it is Hermitian. This 
means that it has a real spectrum, and we may talk about its positivity. The 
positivity constraint placed on the dynamical map earlier turns out to be too 
strong [15], [16] . This in itself is a rich topic, and still a highly debated issue 
in the community [T71 [18] . We will not go into the details of that discussion 
here; though it will be fruitful to define positivity classes associated with maps. 
There are three classes associated with the positivity of dynamical maps: 

Completely positive: A map is called completely positive if all of its 
eigenvalues are positive semi-definite. Such maps always map all positive 
states to positive states [191 |2Q]. Thus, the valid domain for a completely 
positive map is the set of all states. Unitary maps are an example of this 
class. 

Positive: A map is called positive if not all of its eigenvalues are positive, 
but it maps all positive states to positive states. The transpose map for a 
qubit has this quality. It has one negative eigenvalue, but the transpose 
of any density matrix is another valid density matrix [2T] . 

Negative: A map is called negative if it maps any positive state to a 
negative state. These maps have at least one negative eigenvalue. Maps 
of this kind are still physically valid for a certain set of states [22]. This 
set is called the compatibility domain [23] . 
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The set of positive and negative maps are often put in the same category known 
as not- completely positive maps. For our purpose, we will assume that the 
dynamical map need not always be completely positive. In fact the negative 
eigenvalues will come in handy later on in chapter [5j 

2.3.5 Choi representation 

Because the dynamical map has a real spectrum, we can write its action 
in terms of its eigenmatrices {C*-™'''} and eigenvalues {Xm} as 

'Brr';ss'pf.'s' = ^ -^m Crr'Vf's' CL'^ ' (2.22) 
m 

where 1 < m < d"^ (because "B is a (P x matrix, it has d"^ eigenmatrices). If 
the eigenvalues are positive then we can absorb them into the eigenmatrices 
= v^C^™^ to get 

■Bp^ = J2 C^'^^P^ C^"'^\ (2.23) 

m 

with 

J2 cM^c^^") = 1. (2.24) 



Eq. 2.23 is know as the Choi representatior^ This representation was first 



pointed out by Sudarshan et al. [3], however Choi made use of this represen- 
tation to study the positivity classes for maps|^ Any map that can be written 

^ The C-matrices are not unique; they can be unitarily rotated, but the dynamical map 
25 is unique. 

^This representation is also known as the operator sum representation, the Stinespring 
form, and the Kraus canonical form. 
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in the Choi form is completely positive fl9[ [20l EH [2l], which is a convenient 
definition for complete positivity. 

The advantage of writing the dynamical map in the Choi representation 
is that now the dynamics is governed by a set of operators in the space of the 
system. Unitary evolution is a special case in this representation 

p§ ^ Up^U with C(i) = U, C("^) = (for all m > 1). 

2.4 Dynamical maps from contraction 

Let us now consider dynamical maps coming from the contraction of 
Hamiltonian evolution. If we know the initial state of system 



can we find the dynamical map? 

To make our examples physical we consider a toy model by taking a 
combined state of the system and the environment, and letting it develop in 
time unitarily. Then we can find the map by considering the evolution of the 
reduced states of the system 




(2.25) 



and the final state 



(2.26) 



At) 



SpS(to) 



(2.27) 



(2.28) 
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where the initial and final states are obtained by contracting the environmental 
degrees of freedom. This is all we are allowed to know in order find the 
dynamical map. We will workout two examples in two different manners to 
obtain the dynamical map in each case. 

2.4.1 Initially product states 

Let us now construct the dynamical map from physical motivations. Let 
us start by looking at the special case of an initially product state (simply 
separable) p^^ — ® p^. We can write the action of the dynamical map in 
terms of matrix indices as follows 

'Brr';ss'p''r's' = J2 Ure;r'yr's'PapU:,;s'/3- (2-29) 
afSe 

Then we can construct the dynamical map right away by removing p^ from 
both sides. In the S form the map is 

Srr';ss' = ^ Ure-r'aPaj^Ul^-s' p- (2.30) 

The final state of the system is given by 

Prs(^) = ^rr';ss'Pf's'(^o) 

= \Y.^re-yaPipUl^s'Apl'Ato). (2-31) 

We can see in the last equation that the dynamical map acts linearly on the 
state of the system, i.e. it does not depend on any of the parameters of the 
initial state of the system. The linearity here is consequence of the hnearity 
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of quantum mechanics. The map only connects the initial state of system to 
the final state of the system. 

Since is a positive matrix, we can take its square root and distribute 
it. This way we can obtain the Choi representation 



^rr' ^ss' ■ 



Pi 



(2.32) 
(2.33) 



e7 



Finally we have 



-Drr';ss'Pr's' " 2-^ rr' Pr's'^ss' ' 



(2.34) 



This last equation is in the same form as Eq. |2.23[ which proves that maps 
arising from initially uncorrelated states are always completely positive. Con- 
versely, any completely positive map can be thought as coming from a unitary 
evolution of the system and the environment in the product form. 



2.4.2 Initially correlated states 

Let us now look at an example of a dynamical map for a total state 
where the system and the environment are initially correlated. Let us work 
with the simplest open system, a state of two qubits. We will treat the first 
qubit as the system and the second as the environment. The density matrix 
for two qubits in general is 

1 



= ^{I (g) I + ajGj ®l + bkl® (Jk + Cjkaj ® (Tfc}, 



(2.35) 
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where aj are the Pauh spin matrices, aj and bk are the Bloch vector compo- 
nents of the system and the environment respectively, and Cjk represent the 
correlations between the system and the environment. 

For our example we set — Cjk — except C23 to simplify the calcula- 
tions. Then we have the following density matrix for our two qubit system 

p^^ito) = ^{I ® I + ajaj ® I + C23(72 ® 0-3}, (2.36) 

This state is separable but it is not a product state. Let us evolve the total 
state with the unitary operator 

U^e-'"^^ Y[ {cos{u;t)I®I-ism{u;t)aj^aj}, (2.37) 

J=l,2,3 

where 

H ^(j^aj^aj. (2.38) 
j 

The evolved state Up^^{to)U^ is 

p^^[t) = ^{l<^I + aj cos^{2ujt) (7j (8) I - C23 cos{2ujt) sm{2ujt) ai<S)l 
+aj sin^(2a;t) I (8) cTj -|- C23 cos(2a;i) sin(2a;t) I (g) (Ji 
+ai cos{2u!t) sin(2a;i) (cr^ (T„ — (7„ (g) cr^) 
+C23 cos^{2ujt) (72 (8) (73 C23 sin^(2a;t) 0-3 (8) (72} , 

where j runs from 1 to 3 and l,m,n are distinct and cyclic. 

We are now in the position to find the dynamical map for the system 
going from to to t. For that we need the states of the system at those times. 
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These are obtained by tracing the environment out of the total state at these 
two times. The initial state is 



p^(to) = Tr£[p^^(to)] = -{I + a,a,}, 



(2.39) 



and the the state of the system at time t is given by 

p^t) = Tr, [Up''{to)U^] 
1 



{I + cos^ {2ujt) a.ja.j — C23 cos {2ujt) sin {2u}t) ai}. (2.40) 



We can see how the linearly independent elements, that span the space 
of the system, are mapped from to to t 



A : 



ai — ia2 
ai + ia2 



/ 1 + C^as \ 

C^(ai - ^2) - C23CS 
C\ai + ma) - C23CS 

V 1 - c^as y 



(2.41) 



where C = cos{2ujt) and 5* = sin(2ti;t). Since the stochastic map A acts on 
the state as a matrix on a vector, we can construct it by inspection as follows 



A 



1 + 










—C23CS 


2C^ 





— CasCS* 


-C23CS 





2C2 


-C23CS 










l + C 



(2.42) 



Notice that A is not Hermitian, and the element of the first and the last rows 
add to unity. 



more general prescription of one qubit system and one qubit environment is given in 
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To construct the dynamical map, we simply rearrange each row of A 
ling blod 
/ 1 + 



into the corresponding block matrix (see Eq. 2.11) 



^1 n 1 - 

S = - 
2 

\ -C2zCS 1 + C" / 



1-C2 -C23CS 
-C2zCS l-C^ 



(2.43) 



Note that B is Hermitian, the 2x2 block diagonal elements add to unity, and 
is an affine transformation [26] that squeezes the Bloch sphere of the qubit 
into a sphere of radius cos^(2ci;t) and shifts its center by C23 cos(2co't) sin(2co't) 
in the ai direction. 

The eigenvalues of the map are 
1 



Al,2 

-^3,4 



[1 — COS i2u}t) ± C23 cos(2u;t) sin(2u;t)]. 



1 + cos2(2cjt) ± cos(2u;t)A/4cos2(2cjt) + 0^3 sin2(2c<;t) 



It is easily seen that A3^4 are always positive. While, for Ai^2 to be positive, 
we need sin^(2ti;t) > ±023 cos(2ti;t) sin(2co't). We can choose C23 such that 
this condition will be violated for some values of ut making the map S not- 



completely positive. The eigenvalues are plotted in Fig. 2.3 



It has been previously shown that not-completely positive maps come from 
initial entanglement [21] • This example shows that initially correlated states 
can lead to not-completely positive maps. A similar example has been worked 
out in [TTIHH]. The map S has physical interpretations as long as it is applied 



25 




Figure 2.3: The eigenvalues of the dynamical map in Eq 2.43 are plotted as 
function of 2u)t. One of the eigenvalue is negative for certain values of ut; we 
have taken C23 = 0.5. The negative eigenvalue is due to the initial correlations 
between the system and the environment. 

to initial states p^(to) that are compatible with the total state p^^(to) [23]. Let 
us discuss these ideas in more detail. 

The state p^(to) depends on the parameters aj. For p^(to) to be pos- 
itive the Bloch vector must be smaller than unity in magnitude, |ap < 1. 
Additionally, the Bloch vector is constrained by the positivity condition of the 
total state p^^(to)- This constraint arises from the presence of the correlation 
term C23. Because of this term, p^(to) cannot be any state, for instance a pure 



state (see Fig. 2.4 for a graphical explanation). 

The value of the correlation can range from — 1 < C23 < 1. If we fix 
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Figure 2.4: [Left] The total state of the system and the environment has 
correlations, therefore the set states that are physically available for the system 
part (green ellipse on bottom) are less than set of all state (dotted blue circle). 
This is the compatibility domain. [Right] Conversely, if one demands to embed 
set of all system states into the total state (compatible are in the green ellipse 
and incompatible are represented by the blue part of the circle on the bottom), 
then the total state may become unphysical (fuzzy blue and red area above). 



C23 in that range, then for the eigenvalues for p^^(to) to be positive we get the 
following condition 

1> |aip + (|a2| + |c23|)' + |a3|'. (2.44) 

According to the inequality, when |c23| > 0, p^(to) cannot be a pure state; 
otherwise p^^ will not be positive. The volume enclosed by this inequality is 
the valid domain of the map, and the action of the map on this set of states has 
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a physical interpretation. This domain is known as the compatibility domain 
P7] . The map acting outside this domain may transform a physical state to a 



non-physical state as shown in Fig 2.5 




B 





P''(t) 



Figure 2.5: [Top] The total state evolves unitarily, which changes the po- 
larization of the system, the environment, and the correlations between them 
(note the color difference). [Bottom] The dynamical map acting on the com- 
patibility domain results in physical states (green ellipses from left to right). 
But its action on set of all states (dotted blue circle on left) can take physical 
states to non-physical states (dotted red ellipse on right). 



There is another domain associated with dynamical maps. This domain 
comes from the unitary operators, and it called the positivity domain. The 
positivity domain is the set of positive states that get mapped to a positive 
state by the dynamical map. The compatibility domain is smaller than the 
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positivity domain, and it is a subset of the positivity domain. We will not 
show an example of the positivity domain here; it may be found in |23j . 

2.5 Size of the environment 



The Choi representation shown in Sec. 2.3.5 tells us that only a finite 
number of operators govern the most general dynamics of the system. Then we 
may ask what is the largest dimension of the environment that is necessary to 
simulate the most general dynamics? It is often assumed that the dimensions 
of the environment is very large compared to the system. We show here that 
this is not the case. 



In Sec. 2.4.1 we showed that any completely positive map can be 



thought of as a contraction of initial product states evolving. Let us re-examine 



Eq. 2.30 



'S>rr';ss' — ^ Ure-r' aPaP^le^s' 0- (2.45) 
ea/3 

Let us consider the state of the environment to be a pure state 

'Brr'.ss' = $^f/r-.;.'a|0)j0|^t/.V;3 (2-46) 

eQ/3 

= {^Ur.,r>o.\'P)}j{Y.Us.,s'p\4>))j (2.47) 

= (2.48) 

If S is defined on a d dimensional system then S is a x matrix. Hence 
it has only (P C-matrices, therefore the index e only needs to run to d?. The 
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index e also represents the dimensionality of the environment. 

Conversely, any completely positive map can be thought of as a con- 
traction of a pure state of the environment of dimension dP and unitary U 
(see pij for more details). The unitary operator can be constructed from the 
C-matrices as 

C,7, = Ure,r'a. (2.49) 

The value of a is fixed by the choice of the state of the environment. Different 
C-matrices are determined by the value of e. All other elements of the unitary 
operator can be arbitrary, as long as they satisfy the unitarity condition for 
U. 

2.6 Discussion 

Let us review the most essential statements made in this chapter. We 
have shown that dynamical maps govern the most general finite time evolution 
of a quantum state. The not-completely positive maps are due to the initial 
correlations between the system and the environment. These maps have a 
physical interpretation as long as their action is restricted to the compatibility 
domain. 

We also showed that the size of the environment need not be very large 
for the most general evolution. It has to be at most (P, where d is the size 
of the system. With this in mind we will justify the validity of the simple 
examples that will be studied through out this dissertation. 
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It should now be clear that the dynamical map formalism is a very 
powerful theoretical tool. One may ask, are dynamical maps experimentally 
realizable? The answer is yes, quantum process tomography is an experimental 
procedure that allows one determine the open evolution of a system. We will 
review several quantum process tomography procedures in the next chapter. 
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Chapter 3 



Quantum process tomography 



There has been significant experimental interest in the dynamics of 
quantum correlations, entanglement, and coherence in the context of quantum 
information theory [H [2H1 EHl EOl EI] • AH of these categories require the study 
of multi-partite quantum systems. Experimentally, for all of these cases, the 
external influences due to the environment are not uncommon. In the last 
chapter we introduced the dynamical map formalism to describe the open 
dynamics of a quantum system. We now look at the experimental analog to 
dynamical maps, which we call process maps. 

Quantum process tomography [321 |33] is the experimental tool that 
determines the open evolution of a system that interacts with the surround- 
ing environment. It is the tool that allows an experimenter to determine the 
unwanted action of a quantum process on the quantum bits going through 
it. It is an important tool for quantum information processing. A state going 
through a quantum gate or a quantum channel will experience some interac- 
tions with the surrounding environment. Quantum process tomography allows 
the experimentalist to distinguish the differences between the ideal process and 

^Sometimes the process is known, but we are only interested in characterizing the effects 
due to an unknown process. 
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the process found experimentally. Therefore it is an important tool in quantum 
control design and battling decoherence (loss of polarization). 

Today there are many variations of the original quantum process to- 
mography procedure, namely ancilla (entanglement) assisted process tomog- 
raphy ISl |35l |36l [37], direct characterization of quantum dynamics [381 139] . 
selective efficient quantum process tomography^\, and symmetrized charac- 
terization of noisy quantum processes^\. Some of these procedures have been 
experimentally tested [121llSlllllllSlllSlllZlllHllin]. 

3.1 Quantum process tomography: general description 

The objective of quantum process tomography is to determine how a 
quantum process acts on different states of the system. In very basic terms, a 
quantum process connects different quantum input states to different output 
states: 

input states process —>■ output states. (3.1) 

The complete behavior of the quantum process is known if the output state 
for any given input state can be predicted. 

A quantum process can be anything from controlled unitary evolution 
to dissipative open evolution (or in most cases some combination of the two) 
experienced by the state of the system. The quantum process, thus can be 
described by a map just like the dynamical map from the last chapter. The 
only difference being, this is now to be done experimentally. 
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The tomography aspect of quantum process tomography is to use a 
finite number of input states, instead of all possible states, to determine the 
quantum process. For instance, to determine the dynamical map we only need 
to know the mapping of the elements of the density matrix from an initial time 
to a final time, 

S : Pij{to) Pij{t). (3.2) 
The elements of the density matrix linearly span the whole state space. 

Experimentally, we do not have the access to the individual elements of 
the density matrix; we can only prepare physical states. Thus, a set of physical 
states that linearly span the state space will be sufficient for the experiment. 
A state space of dimension d requires (P states to span the space. Once the 
evolution of each these input states is known, by linearity the evolution of any 
input state is known (see P for detailed discussion)]^ 

For example, the following four projections as input states are necessary 
to linearly span the whole state space of a qubit: 

p(^'-) = i(I-aO, P(^'+) = ^(I + aO, 

P(2,+) = 1(1 + ^2), P(3'+) = ^(I + as). (3.3) 

Any state of a qubit can be written as a unique linear combination of these 
four projection^ 

^The linearity of the quantum process is an assumption here. 

^The linear combination will not always be convex. For example P^^^^) = p(i^+) _(- 
_ pi^'+\ Also notice that these four states form a linearly independent set, but they 
are not orthogonal to each other. 
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Using the set linearly independent input states P^"^\ and measuring 
the corresponding output states Q*^™]^ the evolution of an arbitrary input 
state can be determined. Let A be the map describing the process, which we 
call process map, and an arbitrary input state be expressed (uniquely) as a 
linear combination 'YljPmP^^^- The action of the map in terms of the matrix 
elements is as follows: 

r's' \ j / j 

3.2 Standard quantum process tomography 

In quantum process tomography, it is often assumed that at the be- 
ginning of the experiment the state of the system and environment are un- 
correlated. At that point they evolve in a closed form under some unitary 
transformation. Under this assumption we write down the dynamical equa- 
tion in terms of matrix indices 



Qif = J2Ur.yaP^:)p^,U:,,s', (3.4) 

e,a,l3 

J2 Ur.yaP^pu:,.,s'A p^r:) (3.5) 

Arr'-ss'Pri^) ■ (3.6) 



We call Eq. 3.6 the linear process equation. An equation that relates the 



input states to output states is called a process equation. Notice that the 



*This requires performing quantum state tomography, see App. |Aj 
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linear process equation looks very much like Eq. 2.31 which describes the 
dynamical map for the initially product states in Sec 2.4.1 

Experimentally, we do not have any information about the global uni- 
tary operators!^ or the state of the environment. But if we know the output 
states corresponding to the necessary input state, then we find the map by the 
following expression. 

W = ^E^-^^^v*, (3.7) 

n 

where d is the size of the system and P^"^ are the duals of the input states 
satisfying the scalar product 

p(m)tp(n) _ p{m)* p(n) _ r 



The duals for the projections in Eq. |3.3| are 

p(2.+) = a2, P(=^'+) = ^3. (3.8) 

This procedure is known as standard process tomography. 

The map in Eq. (3.7) determines the open evolution of any state sent 



through that process. This is a powerful statement, yet this procedure has 
two serious downsides to it. First, the procedure assumes that a set of linearly 
independent states can be prepared for any experimental setup. This may 



^Unitary operators that acts on a bipartite state are called global, while unitary operators 
that acts only on the subpart are called local. 
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not be the case in reality. In fact for certain situation the experimentahst 
may only care about what happens to certain elements of the density matrix 
and may not care about what happens to the rest of the state space. In 
that case, determining the whole process map is an over kill. For one qubit 
this may not seem like a serious issue, because performing additional two or 
three experiments is generally not difficult. But consider an experiment with 
n qubits. In that case, to determine the whole map requires preparing 2^" 
input states. This may not be practical for all experimental setups (see for a 
detailed analysis [SO])- The second issue is precisely the number of required 
inputs. This number grows exponentially, and so do the number of necessary 
experiments with it. There are elegant methods that assist in overcoming 
both of the issues. We tackle the first issue first with ancilla assisted process 
tomography in the next section. 

3.3 Ancilla assisted quantum process tomography 

Let us consider a situation where the experimentalist has an ancillary 
system available in addition to the system of interest. The ancillary system 
can be used to help overcome the limitations of preparing a set of linear in- 
dependent states. In fact, we will only need to prepare one state with this 
procedure, known as ancilla assisted process tomography [361 137]. Let us once 
again assume that the initial state of the system, environment, and the ancilla 
are uncorrelated at the beginning of the experiment. We start by entangling 
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the ancilla with the system with unitary V that acts on the combined state. 

p-^§ = yp^®pV, (3.9) 

where p-^ is the initial ancilla state, is the initial system state, and R-^^ is 
the maximally entangled state of the system and the ancilla. 

At this point we send the system state through the quantum process. 
We then make a measurement J^"^^ on the ancilla after the system has gone 
through the complete quantum process, and analyze the output state of the 
system. The linear process equation now looks as 

Q^"^^ - ^[j(l)^.t8] Tr^ [^^"^Tre [Up^YU^ J^^^] , (3.10) 

where we have normalized the final state by dividing the probability that the 
ancilla part of p"^^ will collapse to J^"^\ Since the system and the ancilla are 
maximally entangled, we are guaranteed the perfect knowledge of the state of 
the system before it goes through the quantum process by simply knowing the 
state of the ancilla. 

Notice that J^™^ commutes with the unitary U and trace over the en- 
vironment since only the system is going through the quantum process. We 
can pull J*^'"-* inside the trace over the environment 

/^(m) \ \ ^ TT j{m) AS T(m) £ tt* 

~ Trfj(™)o-^^l ^ '-'re;r'a'Jzx h>xr';ys"^yz Ha/S^ se;s' 13 

— \^ TT 7'Mp("*)^£ TT* 

— 2-^ 'J re\r' cx'J zz -^r's' Pa/S'-' se;s'l3 

e,a,/3,z 




where P^*") is the state the system collapses to when the ancilla collapses to 
pirn) jg g^jgQ ^j^g desired input state. Since acting with J^'"^^ on the 
ancilla before the system goes through the process or after is the same, we can 
chose the any input state p(™) by making measurement on the ancilla after 
the system has long passed through the process. Notice the process map A 



has the same form as in Eq. 3^, therefore the rest of this procedure is the 
same as the standard process tomography procedure from the last section. 

By this method we have eliminated the problem of having an exper- 
imental setup equipped to prepare a variety of input states. Now we need 
to prepare only one entangled state of the system and the ancilla. The cor- 
relations between the system and the ancilla do the rest of the work for us. 
This is an example of quantum parallelism. It is interesting to note that this 
method works even when the system and ancilla aren't entangled [36j, though 
not as well. There are quantum correlation other than entanglement that can 
be used as resources as well. A measure of correlations due to Ollivier and 
Zurek and Hederson and Vedral called quantum discord [511 ESI [53] gives some 
insight on how entanglement, quantum correlation, and classical correlations 
are different. 

Though, in solving one problem we have added another. On the one 
hand, our ability to prepare a verity of input states has increased dramatically, 
on the other hand we now have to worry about keeping the ancilla completely 
isolated from the environment. This is almost never practical. We will not 
discuss this particular problem in this dissertation (See. VIII in [5J for a 
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discussion). 



3.4 Direct characterization and other methods 

The trick with the ancilla assisted tomography method is to create a 
maximally entangled state to prepare all possible inputs simultaneously. In a 
sense entanglement is being used as a resource in this method. We can further 
utilize this resource to decrease the number of experiments necessary to find 
the process map using error correction techniques. This is known as the direct 
characterization of quantum dynamics [3HI ES] • 

The basic idea is to modify the ancilla assisted procedure by making a 
joint measurement on the system and the ancilla after the system has passes 
through the quantum process, instead of making a measurement on just the 
ancilla. 



Amn — Trg 



above p-^^^""^ an entangled state (not always maximal) of the system and an- 
cilla. J^™-* now represents a joint measurement on the combined state of the 
ancilla and the system. It turns out that what arises is simply an element 
of the process map A. By choosing proper inputs and proper measurements 
we can find the map without every looking at the output state. The philoso- 
phy here is that determining the process map is important not the input and 
the output states. Therefore in this method quantum state tomography is 
unnecessary. 
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This method gets its inspiration from quantum error correction codes, 
which is beyond the scope of this dissertation. We will not go into the details 
of any of these methods, since they are of no consequence to our work here. 
Though we would like to point out that these methods resolve the issue of 
preparing a large number input states. In the direct characterization of quan- 
tum dynamics the number of inputs necessary to determine the process map 
grow polynomially with the size of the system rather than exponentially as in 
the standard quantum process tomography procedure. This makes it possible 
to carry out quantum process tomography for large systems. There two other 
techniques along this line, selective efficient quantum process tomography |40j . 
and symmetrized characterization of noisy quantum processes pij, that we will 
not review here. 

3.5 Discussion 

In every quantum process tomography procedure above the input states 
are thought to be pure states. There are two advantages of using pure states 
as inputs. First, it is easier to span the space of the system with a set of 
pure state than it is with a set of mixed states. Second, pure states are always 
uncorrelated, which is one of the central assumption in every procedure above. 
This is often called the weak coupling assumption. It is simply a matter of 
preparing necessary pure states to perform a quantum process tomography 
experiment. 

In reality, the weak coupling assumption may not be true for some 
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experiments. In that case the initial state of the system will not be pure. 
Thus, some preparation procedures must be applied on the system to prepare 
a pure state. We relax the weak coupling assumption in the next chapter and 
analyze state preparation for open systems. In chapter |5] we will return to 
quantum process tomography armed with the knowledge of how preparation 
procedures work. 
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Chapter 4 



Preparation of input states 



In this chapter we study the effect of the preparation of the input states 
for a generic quantum experiment. In quantum process tomography, a hnearly 
independent set of states that span the system space have to be prepared. But 
for a generic experiment, we may need to prepare a set that is larger than the 
hnearly independent set. In this section, we will not concern ourselves with 
quantum process tomography and keep the discussion focused on preparation 
only. 

We can describe a quantum experiment in three general steps. The 
experiment begins with an unknown state that has to be altered into known 
input state. After the preparation, the prepared state is subjected to some 
quantum operation, and finally the outcome is analyzed. In this chapter we 
will only concern ourselves with the first step of state preparation. 

Preparation procedures are very complicated in practice. Since we can- 
not describe each preparation procedure in detail, we will attempt to develop 
a general theory of preparation. For clarity, we will first deal with the prepa- 
ration of a general state p, and then embark onto the preparation of a state of 
a system that may be correlated and may interact with an environment. 
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4.1 Preparations in quantum mechanics 

We can summarize a set of preparations by considering an experiment 
that starts with a generic state of the system p, which is then ahered into a set 
of desired inputs {P*^™^}. We now have to somehow connect the initial state 
to the input states. 

As we saw in chapter [2| the most general dynamics of a quantum state 
are described by a stochastic map. We can denote the procedure for preparing 
the mth state with a map [p(™). The only restriction we put on the preparation 
map is that it be completely positive. This is because no matter what the 
initial state of the system is, we should be able to prepare the state. If 7 
is completely positive, then its action is defined on all possible states of the 
system. 

The action of 9^"^-' onto p can be written as 

'p('^) n , , - pi"") (A 1 ^ 

^rr';ss'rr's' — ^rs ■ \^-^) 

Before we look at the preparation procedure for an open system, we 
should discuss how a most general preparation can be carried out with the 



aid of an ancillary system. As we saw in Sec. |2.5[ to describe the most 
general completely positive dynamics of a quantum state of dimension d, we 
only need an environment of dimension (P. Our goal here is to be able to 
prepare the system, so we need control over the system. Therefore, instead of 
an environment we place two ancillary systems of dimension d to perform the 
most general input state. 
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4.1.1 Trace preservation 

It turns out that not all preparation procedures are trace preserving. 
The action of such maps is written as 



r 



J_rnM n , , - PM (A 2\ 

,(m) rr'^ss'h'r's' — ^rs ' l^-^J 

where the normalizing factor 

_ Tr[T(™V] (4.3) 

is the probability with which p will become 

However, there is a complete map that preserves the trace. Consider a 

completely positive preparation map, J", acting on p. We can rewrite this map 



in the Choi representation (Sec. 2.3.5) as 

y(p) = ^7^(™Vv^('")^ (4.4) 

m 

where tt*^™^ are the eigenmatrices of CP satisfying = I, which 

guarantees trace preservation. We now combine, for each m, the eigenmatrices 
7j-(»T*) -with its Hermitian adjoint tt*^™-'^ to get 

y(p) = y^'") p. (4.5) 

m 

It is now clear that the mth preparation is given by the action of y^™), which 
does not preserve trace. But a set of preparations belonging to T, to- 

gether preserve the trace. 
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This suggests that, to prepare m input states, we will in general need m 
preparation maps. In other words, an experiment will require m preparation 
procedures. Not all of these procedures will be completely unrelated, in fact, 
most will share some common features. 



4.2 Preparations in open quantum mechanics 

How does the situation above change when we consider open quantum 



systems? We simply follow the procedure laid out in section 4.1, but instead 
of using the generic initial state p, we must use a bipartite state of the system 
and the environment p^^. The preparation map still only acts on the state of 
the system. 

The preparation of a system belonging to a bipartite state of the system 
and the environment we get 

g£{m) _ _J_^{m) „ S£ (A 

^r'a;s'f3 ~ ^(m) r'r" -,8' s" aa;f3l3Pr"a;s" 13^ 

where indices a and (3 belong to the state of the environment and R^^^""^ is the 
bipartite state of the system and environment after the preparation. Notice 
that the preparation map still acts on only the system part of the initial 
bipartite state. For completeness, in the equation above we have included an 
identity map 3 acting on the state of the environment, but for simplicity we 



will omit writing it from here on. Fig. |4.1| shows graphically the action of the 
preparation procedure. 

Even though the preparation procedure only acts on the state of the 
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Figure 4.1: At the beginning of the experiment an unknown state of the system 
and the environment is present. The mth state of the system is then prepared, 
which diminishes the correlations between the system and the environment. 

system, the state of the environment in general will indirectly be affected by the 
procedure. Suppose our goal is to prepare an uncorrelated state. Simplifying 
the notation above we get 

^S£,(m) _ rp(m)^S£ _ p{m) ^ p£.(™)_ {'^■7) 

In the last equation, we see that the state of the environment has picked 
a superscript m, because if p^^ is initially correlated, then an action on the 
system part will necessarily effect the state of the environment. However, if 
the weak coupling assumption is retained, then the preparation procedures for 
the closed and the open cases are the same. Let us investigate this through two 
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specific procedures. See [51] for a similar analysis for quantum measurements. 

4.3 Two common preparation procedures 

Uncorrelated states are often desired for experiments. Therefore we 
will only focus on preparing pure states. We describe two of the most common 
preparation procedures to prepare pure states (thus uncorrelated) used today. 
We first analyze each preparation procedure for closed quantum systems, and 
then analyze them for open quantum systems [5]. In this manner, the dif- 
ferences between two will be apparent. More complicated procedures can be 
easily studied by writing down the corresponding preparation map. 

4.3.1 Stochastic preparation 

Many quantum experiments begin by initializing the system to a spe- 
cific state. For instance, in the simplest case, the system can be prepared to the 
ground state by cooling it to near absolute zero temperature [HI HHl SZl SHI SH] • 
Mathematically, these set of operations can be written as a pin map [55] . Con- 
sider the following map 

err-ss' = mmrs^r's'. (4.8) 

The action of this map on the initial state of the system is 

r's' 

= 1$) ($|Tr[p] = 1$) ($1 . (4.9) 
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In the procedure above, no matter what the initial state of the system was, it 
is "pinned" to the final state |$) ($|. Along with the pin map, a set of unitary 
transformations on the system are sufficient to prepare the mth input state. 
In that case wc can write the preparation as a composition of the stochastic 
pin map and a set of local unitary transformations (acting only on the system) 



We call this procedure the stochastic preparation method. 

Each preparation map is trace preserving, because the pin map and the 
unitary transformations that follow it are all trace preserving transformations. 
Also, to prepare m input states, we need m unitary rotations, hence we have 
m preparation procedures. However as claimed before, each of these procedure 
has the pin map as a common step. 

For a closed system, we could just as well use m pin maps to prepare 
m input states. The trouble only arises when the system is interacting with 
an environment. Let us look at this in detail. 

4.3.1.1 For open systems 

Let us consider the action of the pin map © on a bipartite state of the 
system and the environment. 



as 



(4.10) 



1$) ($|®Trs[p§^] = 1$) 



(4.11) 
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As we can see, not much happens differently when a stochastic map is apphed 
to a bipartite state. The only difference is the presence of the environmental 
partner . The pin map fixes the system into a single pure state, which 
means that the state of the environment is fixed into a single state as well. 
The purpose of the pin map is to decouple the system from the environment, 
to eliminate any correlation between the system state and the environment 
state. 

Once the pin map, ©, is applied, the system is prepared in the various 
different input states by applying local (only on the system part) unitary 
transformations as before. The prepared state Bf'^ takes the following form. 

Strictly speaking, in a real experiment the pin map is a complicated 
set of steps. For instance, cooling the system so it relaxes to the ground state 
is very different from optically pumping a system to an excited state. Two 
different pin maps would correspond to these two different methods. Then, the 
state of the environment in each case can also be different. We have assumed 
that an identity map acts on the environment in our mathematical formulation. 
This may not be the case; we cannot analyze the state of environment before 
or after the pin map is applied. Therefore, to be careful, we should label as 
p% to clarify that the state of the environment may depend on certain steps 
involved in applying the pin map. 
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Then, it should be emphasized that the initial pin map is critical. It 
may be tempting to simply use a set of pin maps, e^™), to prepare the various 
input states p('"). The different input states will be 

^S^C") ^ rp{m) SE ^ @(m) S£ 



As seen in the last line of the equation above, the state of the environment 
has taken the subscript B''™^ through different pin maps along with the input 
state. What this means is that the state of the environment is affected by 
the preparation procedure. We show in the next chapter, when the prepara- 
tion procedure is not carried out carefully in a quantum process tomography 
experiment, the process map can behave non-linear. 

4.3.2 Projective preparation 

Let us look at another preparation method that is also very common. 
Consider an optical experiments, where a beam of photons is sent through a 
polarizer |l2l SSI HI]. Polarizers project the polarization vector of the photons 
into a set orthogonal directions. By using a quarter wave plate, all necessary 
input states can be prepared. 

We can write the action of mth preparation map acting on the initial 
state of the system as the following: 



(4.12) 



yMp 1 



pi'm.) pp{m)^ — p{m) 



(4.13) 



Tr[y(™)p] rM 
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where ^^^^ is the normahzation factor with A™'^ = Tr[T*^'"^p] = Ti[P^™'^p]; this 
is the probabihty of obtaining that particular input state from a von-Neumann 
measurement Since projective transformation are not trace preserving, 

the preparation map does now preserve trace. The normahzation factor, ^^j^, 
is simply the inverse of the probability that the system will collapse to the 
mth state, which can be measured experimentally. In this fashion we can 
prepare any pure state as long as we properly normalize them. We call this 
method projective preparation. Note that to prepare m input states m different 
preparation procedures are necessary. 

Though each projective preparation map does not preserve trace, there 
is an overall map that preserves trace. Let us write a complete map in the 
operator sum representation, where we let the eigenmatrices of the map to be 
a set of orthogonal projections. Thus we have 



where {p(™)} form an orthogonal set satisfying p(™-)p(") = 5„„p(™-) and 
j2^pim)^p{m) ^ ^^pM = I. Now we can define each map for a 
single preparation as p(™')p(™)^. Now we see that for a non-trace preserving 
preparation there is a complete map that does preserve the trace. In that 
sense, the projections are the eigenmatrices of the map 

4.3.2.1 For open systems 

A projective preparation is given by a set of orthonormal projections. 
If a particular outcome is observed from the measurement, the state of the 




(4.14) 



m 
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system collapses to that corresponding projection. Therefore, the input states 
can be prepared for an experiment by suitably fixing our measurement basis. 
With the knowledge of the basis and the outcomes, the exact input state 
becomes known. 

Suppose the mth input state is prepared, given by the projection p("*), 
by projecting the system into that state. For the open system, the total state 
becomes 

where once again r^"*^ = Tr[CP^"^V^^]Tr[-P^"^V^^] is the normalization factor. 
Notice that the state of the environment has picked the superscript m. This 
is due to the fact that in general p^^ contains correlations between the system 
and the environment. We should once again emphasize that the preparation 
procedure only acts on the system part and affects the state of the environment 
only indirectly. 

The dependance that the state of the environment picks up is very 
easy to demonstrate. Consider the following one qubit system and one qubit 
environment state: 

p§£ = !(][ (g) I + ai (g) (71 - (T2 ® 0-2 + 0-3 (g) (73). 

This is a maximally entangled state of two qubits. It is easy to check that the 
projective preparation of this state with different projections p(='=-?) — |(I±(Tj) 
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for {j = 1, 2, 3} leads to different states of tlie environment in each case: 



1 



_p(±i)^S£ p{±j) 



r(±i) ' 
p{±j) ^ pi±j)_ 



(4.16) 



Once again, we see that the state of the environment is affected by the prepa- 
ration of system state. This dependance cannot be ignored, and should be 
accounted for whenever the system evolves openly. In the next chapter, we 
will show that when the preparation procedure is ignored in quantum process 
tomography, the process maps obtained can behave non-linearly. 




Figure 4.2: The figure above shows how that for two different projective 
preparations lead to the desired input states, but the corresponding states of 
the environment can be different. In both cases the initial total state is the 
same. 
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4.4 Discussion 

The examples above only aspire to prepare pure states. However, much 
more can be obtained from the prescription spelled out above. In principle, 
nothing is preventing us from preparing mixed states by using so called gen- 
eralized preparations (see [3]). Furthermore, our discussion so far has largely 
focused on preparing uncorrelated states; the very same techniques can be 
applied in preparing entangled states as well. 

We should also add that any time an ancillary system is used in an open 
experiment, it is an unreasonable assumption that the ancillary is uncorrelated 
with the environment. Furthermore, the ancillary may even interact with the 
environment. We will not discuss this topic in this dissertation, but see Sec. 
VIII in [5j for a discussion. 

Note that, if the system evolves in a closed form then all preparation 
procedures are equivalent. Trouble only arises when the preparation procedure 
indirectly affects the state of the environment which then interacts with the 
system. Additionally, even though we will only use the results here in ana- 
lyzing quantum process tomography procedures, they apply to any quantum 
experiment that interacts with an environment. 
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Chapter 5 



Preparation and quantum process tomography 



Each of the quantum process tomography procedure discussed in Chap. 
[3] assumes the system and the environment are initially uncorrelated. What 
if we depart from this assumption? In chapter [2] we showed that initially 
correlated states can lead to not-completely positive dynamics for the system. 
In many recent experiments, the process maps that characterize the quantum 
operations have been plagued with negative eigenvalues and occasional non- 
linear behavior (see references in [5]). We now examine the how the two cases 
are related. 

The quantum process tomography procedures we reviewed require the 
input states be uncorrelated with the environment. Just before the experiment 
begins, in general the state of the system may be correlated with the environ- 
ment. Thus, at the beginning of the experiment it is necessary to prepare 
the initially correlated total state into an uncorrelated state. The difference 
between the process maps found from a quantum process tomography exper- 
iment and the dynamical maps calculated theoretically is precisely the act of 
preparation of input states. Let us investigate this issue by analyzing the steps 
involved in a quantum process tomography experiment. 
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5.1 Quantum process tomography experiment in steps 

The basic steps in a quantum process tomography experiment are bro- 
ken down below: 

• Just before the experiment begins, the system and environment is in 
an unknown state p^^. The system and the environment in general are 
correlated; we are no longer making the weak coupling assumption. 

• The system is altered to a known input state by a preparation procedure. 
The system and environment state after preparation is therefore given 



The input state is given by taking trace with respect to environment 



• The system is then sent through a quantum process. We consider the 
evolution to be a global unitary transformation in the space of the system 
and the environment: 



• Finally the output state is observed. Mathematically it is the trace with 
respect to the environment 



by 



1 



rp(m) S£_ 




1 




(5.1) 



we will call this the generalized process equation. 
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• Finally using the input and the output states, we construct a map de- 
scribing the process is constructed. 

The procedure above is identical to the procedure to find a dynamical map, 
except for the preparation procedure. In the next two sections we will analyze 
quantum process tomography with the two preparation procedures discussed 
in last chapter. The differences between what we find here and dynamical 
maps will be due to the preparation procedures. 



5.2 Stochastic preparation 

In this section we analyze quantum process tomography procedures us- 
ing the stochastic preparation procedure to generate input states. The stochas- 
tic preparation map's action on a bipartite state is given by |4.11 We plug 



that into the generalized process equation 5]T]to get 



= Tr£ [f/p("^) ® p^t/t] . (5.2) 



The last equation is the same as the linear process equation (Eq. 3.6) from 
Chap. [3j Therefore once the input states are prepared, the procedure for 
quantum process tomography is the same as the given in Chap. |3} This is a 
generalization of linear quantum process tomography. Here, we have assumed 
that only one pin map, 6, is used. Therefore we have not labeled the state 
of the environment with the pin map. We will relax this assumption in the 
second example below. 
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5.2.1 An example with stochastic preparation 

Consider the following two qubit state as the available state to the 
experimenter at t = 0_: 

= 1 (I ® I + ajaj O I + 023^2 ® (Tg) . (5.3) 

This is the same state we used in Chap. [2] for the dynamical map example. 
Let the pin map, G, for our example be 

6 = 10) (01® I, (5.4) 

where |0) (0| is a pure state of the system. The preparation of p^^ with this 
pin map leads to 

epS^ = 10) (01 ® h, (5.5) 

yielding the initial state |0) (0| for the system qubit and a completely mixed 
state for the environment qubit. The next step is to create the rest of the 
input states using maps In this case, the fixed state |0) (0| can be 



locally rotated to get the desired input state p(™) (given in Eq. 3.3) 



fi^"^^ 10) (01 ® = V^^'") 10) (01 v^("^)^ ® 

= P(")®^I, (5.6) 

where m = {(1, — ), (1, +), (2, +), (3, +)} and are local unitary operators 
acting on the space of the system. 
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Now each input state is sent through the quantum process. The output 



states can be calculated using Eq. 5.2 Once again let us choose the same 



unitary transformations as in Eq. 2.37 The output states are 



{I - cos^(2cjt)cTi}, Q 

{I + COS^(2cjt)(T2}, Q 



(1,+) 



(3,- 



{I + cos^(2cjt)o-i} 

{I + cos2(2u;t)a3}. (5.7) 



The linear process map is constructed using Eq. |3.7[ the duals in Eq. 

3.8 and the output states: 

/ 1 + C2 2C2 \ 





2C2 1 + ^2 



A. 



1 



(5. 



where C = cos{2ujt). The eigenvalues of the process map are plotted in Fig. 



5.2.2 An example with multiple stochastic preparations 

The pin map used in stochastic preparation must be used consistently. 
Let us show an example of what happens when two stochastic preparation 
procedures are used. Consider a quantum process tomography experiment 
where the following linearly independent states are used to span the space of 
the system 

P' = h, P(l+), P(2+), p(3+)_ (5_g) 
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Figure 5.1: The eigenvalues of the process map in Eq. 5^ are plotted as func- 
tion of 2u}t. As expected the eigenvalues are always positive. The stochastic 
preparation procedure allows the experimentalist to prepare any pure state 
for the system. By convexity, then all possible states of the system are in the 
compatibility domain of the process map in Eq. |5.8[ 



These states form the linearly independent set {I, cTj}; which is different from 



the set given in Eq. 3.3 



Once again, we take the same initial state and unitary operators as in 



the last example. Now suppose the pin map given in Eq. |5.4| is used to prepare 
the state |0) and then by local transformations P^^'^\ p(2"+), and p('^'+) 
are prepared. Finally the mixed state is prepared by letting p('^+) decohere. 



This is a different pin map than the one in Eq. 5.4, so we are using two pin 
maps to prepare two sets of input states. The unitary operator, we are using, 
is often called the swap gate, because it swaps the states of two qubits with 
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the period of t = Then at t = the total state will be 

p^^ (£^) = ^{I®I + a,I®a,}. (5.10) 

The state of the system has fully decohered. 

The corresponding output states for the input states above are found 



using Eq. 5.2, The state of environment in that equation for input = p 
is = ^{I + djCTj}, while for the other inputs the state of the environment 
= |l. Then the corresponding output states are 

Q(') = ^{I + sin'(2cut)a3}, Q^^'^^ = ^{I + cos\2u;t)ai}, (5.11) 
= i{I + cos\2cut)a2}, Q(=^'+) = ^{I + cos\2u;t)a3}. (5.12) 

The last three output states are the same as in the last example, but the fourth 
one is different. 

Now suppose we calculate the output state corresponding to the input 
p(-i). By linearity we have P^"^) = I - P(+^). If the process is linear, then 
the output state for this input state is given by 

g(-i) = 2gm-Q(+i) (5.13) 
= ^{1 + 2sin2(2cjt)o-3 - cos^(2u;t)o-i}. (5.14) 

This state is not physical for certain times, and therefore the linearity and the 
positivity of the process is violated. 

The process has not changed from the last example, only the method 
of determining the process has. Thus when the stochastic map is not used 
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consistently, the process map can behave nonhnear. For completeness we find 
the process map for this example, 


2C2 



A. 






1-^2 



2C2 




\ 



\ 



(5.15) 



where C = cos(2ti;t) and S = sin(2co't). The eigenvalues of this process map are 



negative, as seen in Fig. 5.2, The fact that the state of the environment is not 
a constant for all of the input states leads to the negativity and non-linearity 
of the process map. 

The origin for the negative eigenvalues here is completely different than 



the negative eigenvalues found for the dynamical maps in Sec. 2.4.2 The neg- 
ative eigenvalues in the dynamical map were due to the initial correlations 
between the system and the environment. Since none of the inputs were corre- 
lated with the environment here, negativity arises from the inconsistencies in 

is 



the preparation procedure. Furthermore, the dynamical map in Sec. 2.4.2 



linear and has a physical interpretation within the compatibility domain. No 
such interpretation is found for the process map in this example. 



While this example may not seem realistic, the point regarding inconsis- 
tency arising from multiple pin maps still stands. In realistic cases, the trouble 
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2a)t 



Figure 5.2: The eigenvalues of the process map in Eq. |5.15| plotted as function 
of 2ujt. One of the eigenvalue is negative for certain values of ut. The negativity 
is due to the inconsistency in the preparation procedure, and not due to the 
initial correlations with the environment. 



may not be seen so easily, due to complicated interactions with the environ- 
ment. We will discuss an experiment where multiple stochastic preparation 
procedures are implemented in Chap. |8| 



5.3 Projective preparation 

Let us now consider projective preparation procedure in quantum pro- 
cess tomography experiments. Combining the projective preparations given 
by Eq. |4.15| with the generalized process equation 5J^, we get the projective 
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process equation for the mth input state 



Q(™) = ^Tr£ [f/P^^V^^^^^^f^l • (5.16) 



Is this process given by a hnear map? Dynamically, the evolution of the 
total state p^^ is linear, because every operation on p^^ is linear. However, for 
the purpose of tomography, the dynamics of the prepared input states p(™) 
is of interest. But p(™) appears twice in the process equation, therefore the 
output states, Q^'^\ depend bilinearly on P^*"). This bilinearity can also be 
seen from the dependence of the environment state on the state of the system. 
To see this expand Eq. |5.16 



as: 



gM = TrB [f/P(") ® p^'^'^^t/t] , (5.17) 
with (5.18) 
p^'M = -l^Trs [P(™)p§^] . (5.19) 

The last equation shows the environment state is in effect a function of P^^\ 
It is well known [T71 [STJ [58] that if the initial state of the system is related to 
the state of the environment by some function / as 

pS^ = ^p|®/(pS), (5.20) 

3 

then the evolution of the reduced matrices, p^, cannot be consistently described 
by a single linear map. Above /(p^) are the density matrices of the state of 
the environment with dependence on the state of the system. In this case, the 
function / is of a specific form that gives us a bilinear dependence. 
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It is instructive to now look at a simple example to demonstrate that 
when input states are prepared by projective preparations, the results cannot 
be consistently described by a linear map. 



5.3.1 An example with projective preparation 

Suppose we are unaware of the dependence of the state of the environ- 
ment on the state of the system due to projective preparation and assume that 
the process is given by a linear map. We would prepare a set of linearly inde- 
pendent input states, then construct the linear process map from the output 
states and the duals of input states. 



We once again start with the total state given by Eq. 2.36 and the 



input states given in Eq. 3.3 The state of the system plus the environment 



after each projective preparation takes the following form: 

p(m)^S£pM ^ p{m) ^ Ij (f^^ ^ _ 1^^^ ^^Y^ 

p(2,+)^S£p(2,+) ^ p{2, + ) ^ 1 (I + ^gCTs) , (5.21) 

where Coq = t^^. 

Following the same recipe as before, we obtain the output states using 



Eq. 5.16 We will use the same global unitary from the last two examples 



given by Eq. |2.37 The output states are as follows: 



g(^'-) = ^ {l-cos2(2cjt)ai} , (5.22) 
g(i'+) = ^ {I + cos\2ujt)ai} , (5.23) 
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Q 



(2,+) 



{I — cos(2co't) sin(2co'i(:)cri 



(5.24) 



Q 



(3,+) 



+ cos^(2co't)(T2 + C23 sin^(2u;)f:)(T3} 
{l + sm2(2cjt)(T3} . 



(5.25) 



The linear process map is constructed using Eq. |3.7[ the duals in Eq. 
3.8 and the output states: 



A. 



1 





2C2 + ic'^^CS 







2C2 - ic'23CS \ 



1 - C2 icLCS 



-23'- 

ic'2352 



(5.26) 



Now consider the action of Ap on state p(2'~) = i(l — (T2). p(^'~) = 
p(i.-) _|_ p(i!+) _ p{2,+) jg linear combination of three of the input states used 
above. If the action of A^ is linear then the output state corresponding to 
P(2.-) should be 

A,(P(2H) = Ap (p(i'-) + p(i,+) _ p(2,+)) 

= g(i'-) + g(i.+) - g(2.+) (5.27) 

= ^ (I + c'23 cos(2cjt) sm{2ujt)ai (5.28) 
— cos^(2co't)cr2 — C23 sin^(2ct;t)(T3) . 

Let us check if this the same as if the system was prepared in the state p(^'~) 
by projective preparation. The output state for input P^^-^) can be calculated 
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2a>t 



Figure 5.3: The eigenvalues of the process map in Eq. 5.26 are plotted as 
function of 2ujt. One of the eigenvalue is negative for certain values of ut; we 
have taken C23 = 0.5. The negativity here is due to the initial correlations 
and the inconsistence preparation procedure. 



using Eq. 5.16 



Q 



(2,-) 



- (I — C23 cos(2c(jt) sm{2ut)ai 

— cos^(2c<jt)(T2 — C23 sm'^{2ujt)a3) 



(5.29) 
(5.30) 



where C23 



l-a2 ' 



The output state predicted by the linear process map in Eq. 5.27 is not 



the same as the real output state calculated dynamically in Eq. 5.29, hence 



the linear process map Am does not describe the process correctly. This is not 



surprising; observe the state of the environment in Eq. 5.21 It depends on 
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02, and subsequently the linear process map depends on 02, thus on the initial 
state of the system. The dependence of the environment on the system is the 
cause for the non-linear behavior of the process. 



5.4 Discussion 



The examples with multiple stochastic preparations (in Sec. 5.2.2) 



and projective preparations (in Sec. 5.3.1) look very similar. For multiple 



stochastic preparations, the state of the environment depends on the stochastic 
map. If we think of each projective preparations as an independent stochastic 



preparation then the state of the environment, in the example in Sec. 5.3.1 
depends on the preparation map. Then are the two situations the same? Let 
us show below that the two situations are fundamentally different. 

In the weak coupling limit, the projective preparation procedure does 
not play an important role. As we saw in the last chapter, when the sys- 
tem and the environment are initially uncorrelated, the projective preparation 
procedure will have no affect on the state of the environment. This is be- 
cause the state of the environment is affected only indirectly due to the initial 
correlations between the system and the environment. Thus the projective 
preparations, in the weak coupling limit, yield a linear process map. 

This is not the case in the multiple stochastic preparation example. 
When multiple stochastic maps are used to prepare different input states, the 
inconsistencies do not stem from the initial correlations between the system 
and the environment. In fact, in our example all input states are initially 
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uncorrelated from the environment. The inconsistencies arise from the prepa- 
ration procedures themselves, leading to different states of the environment 
for different input states. These inconsistencies will be absent in the case 
where the system develops in a closed form, since the system does not feel the 
presence of the environment during the quantum process. Hence in the "weak 
interaction limit" multiple stochastic preparations yield a linear process map. 

Lastly, note that each of the process map found in this chapter, given by 
Eqs. 5.8| 5.15, and |5.27| are different from each other. This clearly shows that 
the preparation procedures play a non-trivial role in open quantum system 
experiments. The preparation procedure is the only thing that distinguished 
each case. For the case where no preparation procedure is applied, i.e. the case 



of the dynamical map (Eq. 2.43), the situation is still different. The dynamical 
map, which has negative eigenvalues, is linear and has a valid interpretation 



within the compatibility domain. The process maps in Eqs. |5.15| and |5.27| do 
not have any consistent interpretation. 

When negative eigenvalues are found in a process map, they hint to 
some problem in the preparation procedure. Although this is a bit premature 
to state at this point; we analyze the negative eigenvalues in process map in 
more detail in the next chapter. In Chap. [7] we show that a processes can be 
correctly described using a general dynamical map, which we call dynamical 
M-map. This map is independent of the preparation procedure. 
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Chapter 6 



More on stochastic preparations 



Dynamical maps [21 H] are a generalization of the unitary transforma- 
tion much like density matrices are a generalization of the pure state rays. 
Dynamical maps allow the description of stochastic processes as well as the 
evolution of open systems. Dynamical maps used to describe the evolution of 
open systems are usually defined with a constant environment state p^: 



Implicitly, the state of the environment, p^, is a parameter of the map B. 
Therefore, the linear dynamical map S would only consistently describe an 
experiment if different input states can be prepared independently of the 
environment state p^. The actual issue of how this can be executed is not 
addressed. 

It is possible to consider a dynamical map where the environment is not 
fixed, such as the reduced dynamical evolution of a initially correlated state 



The dynamical map in this problem is applicable only over a compatibility 
domain of states, rather than over the complete state space of the system (see 




(6.1) 



[27]: 




(6.2) 
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Sec. 2.4.2). The compatibility domain is the set of states of the system that 
are compatible with the correlations in p^^. Formally, this problem defines an 
embedding mapQ that relates the initial state of the system to the tota initial 
state, p^^. The embedding map is linear but not necessarily a completely posi- 
tive map [151 [13 EH] • How do we handle the initial correlations experimentally 
in the case of stochastic preparation? 

The stochastic preparation method provides a way for an experiment 
to be performed so that different input states can be prepared with a fixed 



environment state. Consider the process equation 5.2 in section 5.2 



Q{m) ^ r^^^ o Qp^^U^] 

= Tr§ [t/p(™) ® p^'^U^ . 

Then the action of the process map on an arbitrary state of the system is then 
given by: 

ApS = Tr£ [Up^ ® p^'®f/t] . 

Therefore, in this context, the dynamical map is equivalent to the process 
map. However, for consistency, we have to remember that the environment 
state is a constant to the problem, therefore the pin map G should also be a 
constant to the problem. 



^Thc embedding map is also known as extension map, assignment map, or a preparation 
(not to be confused with the preparation map from Chap. |4]). 
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6.1 Pseudo-pure states 

Demanding pure input states in a quantum process tomography exper- 
iment guarantees initially uncorrelated state of the system. Though if it is not 
possible to prepare pure states, we can still determine the process. 

Suppose the initial state is prepared by a pin map leading to 

GpS^ = (pP(3,+) ^ (1 _ p)p(3,-)) ^ J ^Q_3) 
= ^{1 + pas} 01, (6.4) 

where << p < 1. The rest of the input states are prepared by rotations. 
Then the input states are 

= 1{I- pai} P(^'+) = ^{I + pa,} (6.5) 
P(''+) = ^{I + per,} P^''^^ =^-{I + pas}. (6.6) 



For the unitary operator in Eq. 2.37 the corresponding output states will be 



Q^^''^ = ^{I-pcos\2Ljt)ai}, Q^^'+^ = ^{I + pcos\2Ljt)ai} (6.7) 

Q(2.+) = i{i + pcos2(2cjt)(T2}, Q(3'+) = -{I + pcos2(2cjt)(T3}. (6.8) 
2 2 



The only change that we have to make to find the process map is to 
define a dual proper set, in this case 

P^''-^ = ^{pl -ai-a2- as}, P^''+^ = ^{pl + ai - a, - a^}, 

P(2'+) = -a,, P(^'+) = -as. (6.9) 
p p 
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We can find the process map using Eq. 3.7 The process map here case turns 



out to be the same as in Eq. 5.8[ This is to expected; the process map for 



stochastic preparation is defined over set of all states, including the input 
states above. 

The downside of course is that there is no way to distinguish the states 



in Eq. 6^ from the following state 



Gp^^ = -{I^I+pas + CasaaOas}. (6.10) 



Unlike in Eq. |6.5[ this is a correlated state. Yet, both total states have the 
same reduced state for the system part. 

Even in this case we can find the process map properly if we use the 



correct dual set given in Eq. 6.9 The process map in that case will be the 



same as the dynamical map given in Sec. 2.4.2 by Eq. 2.43 In both of the 



examples above, had we assumed that the input states were close enough to 



the pure states we desired, and used the dual set give by Eq. |3.8[ then the 
process map would contain an error of magnitude 1 — p. 

These two examples illustrate the contrary argument to what some 
have suggested, that only the completely positive process maps should be 
considered physically valid. We showed in the second example, that one can 
obtain a not-completely positive process map in a consistent fashion. The 
real issue is that the correlations with environment are not convenient for 
experimental purposes. But that does not mean that one should fix a not- 
completely positive process map to a completely positive process map with 
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numerical methods [601 EI] • 



6.2 Negative maps due to control errors 

There is one more possibihty for non-hnearity and negativity in a pro- 
cess map. This has to do with poor control in the stochastic preparation 
procedure. Suppose the initial state is prepared well using a pin map 

After obtaining this state, the other input states are prepared by local rota- 
tions. Let us consider the case where one of the rotation is not perfect. 

V^'''^ |1) ^ {VT^e\l) - VT+~e |0)) , (6.12) 

where e is taken to be a small positive real number. We introduced a small 
error for the preparation of P^^~\ but we have kept the error simple by not 
giving it an additional phase, i.e. keeping e to be real. For simplicity we have 
also assumed that the error of this sort occurs in the preparation of only one 
state. 

Let us now pretend that we are not aware of this error. Then in reality 
we have the following set of input states 

P^^'") = ^{I + ea,-Vl^a,}, pa.+) = 1{I + aj, (6.13) 
= ^{1 + ^2}, P(=^'+) = ^{I + a3}. (6.14) 



Let us use the same unitary evolution as before given in Eq. 2.37 The output 



states corresponding to the last three input states in Eq. 6.13 are the same as 
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before, given by Eq. 5/I_ For the input state P^^' \ the corresponding output 
state is follows 





^{I + (eas - Vl - eVi) cos^(2cjt)} 




^{I + cos2(2wt)ai} 








I 

-{I + cos^(2cjt)a3}. 



(6.15) 



we can 



Using these output states, the duals given by Eq. |3.8[ and Eq. |3.7 
find the process map. 

The process map turns out to be rather complicated, and its eigenvalues 
are even more complicated looking. Therefore, we do not write them down, 
instead we have plotted the eigenvalues as function of 2ujt. We take the value 



for the error to be e = 0.1 for the plot in Fig 6.1 One of the eigenvalue in Fig. 



6.1 is negative for certain times. This shows yet another cause for negative 
eigenvalues in a process map. The negative eigenvalues here have nothing 
to do with the initial correlations between the system and the environment. 
The negative eigenvalues are attributed to poor control in the preparation 
procedure. 
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Figure 6.1: The eigenvalues of the process map found using the output states 



in Eq. 6.15 and the duals in Eq. 3.8 One of the eigenvalue is negative for 
certain times; we have taken e = 0.1. The negativity is due to the errors in 
the unitary operation implemented to prepare one of the input states. 



6.3 Discussion 



Now we have discussed several scenarios that can lead to negative eigen- 
values for a process map. Two comments are in order at this point. When the 
prepared input states are not pure and the process map has negative eigenval- 
ues, then one should be weary of initial correlations with the environment. To 
further check this, the process map should be tested for linearity by sending 
several additional input states through the process. If the process map pre- 
dicts the output states properly, then one can be confident that the system is 
initially correlated with the environment. If the process map does not predict 
the output states correctly (non-linear behavior), then there are additional 
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problems with the experiment, including the possibility of poor control in the 
preparation procedure. In the case where the input states are pure and the 
process map has negative eigenvalues, the negativity can only come from either 
inconsistencies in the preparation procedure or poor preparation control. 
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Chapter 7 



Dynamical M-map 



As we saw in Chap. [5} the additional step of preparation leads to 
complications for quantum process tomography experiments. At the same 
time, this step cannot be avoided for systems initially correlated with the 
environment. Let us try to attack this problem directly. Let us rewrite the 



process equation 5.1 



qM = -l^Tr£[f/y(™V^^f^1- (7-1) 
In terms of matrix indices we have 



where the sum over e denotes the trace with respect to the environment. We 
are interested in the reduced dynamics of the system states; not so much in 
the details of the preparations. Since the preparation map only acts on the 
system, the trace with the environment has no effect on it. Thus, we can just 
pull the preparation map out of the trace, 

^rs ~ ^(m) ''r'r";s's" / j re;r' cxHr" a\s" (3^ se;s' 13 

e 
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In the last equation, the matrix M is defined as: 



Tyr{'-^) _ 
•^'Vr";s"s' ~ 




S£ 



r"a,s"f3'-^ S€,s'f3- 



(7.4) 



Note that in Eq. |7.3 the superscript indices on M match the elements on the 



left hand side of the equation, while the subscript indices are summed on the 
right hand side of the equation. 

The output state, Q^'^\ is given by the matrix M acting generally on the 
preparation of the input state, Therefore, matrix M fully describes the 

process before any preparation is made. M contains both U and p^^; however 
knowing M is not sufficient to determine U and p^^. As expected, it should not 
be possible to determine U and p^^ through measurements and preparations on 
the system alone without access to the environment. Conversely, M contains 
all information necessary to fully determine the output state for any prepared 
state. We will call M the dynamical 'M-map or just M.-map for short. The 
M-map is due to Aik-Meng Kuah Let us look at some basic properties 
of M-map. 

Trace of JA-map: 

Let us start with the trace of M. To take the trace, we equate indices 
r with s, r' with s', and r" with s" to get: 



Tr[M] = 
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Since WU = ^1^6 ^reypUreya = Ir'a,r'/3, then 

Tr[M]= J2 = ^' 

a,/3r"r' 

where d is the dimension of the system. 

Hermiticity ofM-map: 

As with the case of hnear process maps, matrix M is hermitian. This 
is easy to see by taking the complex conjugate of matrix M, 



S£ TT* 

r"a,s"p'-'se,s'P 



-''Vs";r'V 



The complex conjugate of M is not only the transpose of M, but each element 
of M is also transposed. Hence M is a Hermitian matrix. 

Positivity ofM-map: 

The M-map is composed of a unitary matrix operating on a density 
matrix. Then we can take the square root of the density matrix to get 



•'^W'r",s"s' 



e 

— I Ure,r'a 

€ \ a 



'P 



S£ 



P' 



§£ 



s"|8,(T7 
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where M = U \J p^^ and jji = a'ye. We have written the M-map in operator 

sum representation, hence it is completely positive. This means, the M-map 
acting on any completely positive preparation procedure will lead to a physical 

state. 

The action of M-map can now be written as 



The advantage of dealing with the M-map is that we no longer need 
to worry about the preparation procedure. Once we know the M-map, we 
can act with it on any preparation procedure to find the corresponding output 
state. Since we took out the preparation part out of the M-map, it contains 
all of the dynamical information for the system. In the next section we will 
try to extract some of this information from the M-map. 

7.1 Uninterrupted dynamics 

The M-map contains the dynamics of the system before any preparation 
is made on the system. Perhaps we can learn something about that dynamics. 
Let us start by looking at the definition of the M-map. 



The system is labeled by indices r" and s" . Then by tracing over 
everything else we can find the initial state of the system from M-map. 




(7.5) 




(7.6) 



e 




(7.7) 



rr' 



e 



rr' 
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- ^ Ir'r' ^ lal3plna;s"P (7-9) 
r' a/3 

= cipf„,„. (7.10) 
Furthermore, if we apply an identity preparation, defined as 

on the M-map then we get the final state of the system. 

'yr(rs) fp(I) _ TT Tvf IT* (7 1 9^1 

-''Vr";s'r"'^r'r";s's" ~ %''-r"-^'"-r'r";s's"^s"s' \'-^'^) 

= Ure,r'aK'r"Pf^a-s"0^s"s'Use,s'0 C^'^S) 

= Ure,r'apfux;s' iB^se^s' (7-14) 

= P'.W- (7.15) 

Knowing the initial and the final states of the system (without any 
preparations) allows us to find the dynamics uninterrupted. From this infor- 
mation we can find the mapping from the initial to the final state. In general 
this will not be the same as the dynamical map of Chap. |2} The dynamical 
map is defined over a set of initial states of the system. Here we get a mapping 
from an initial state to a final state. 

7.2 Memory due to correlations 

When a system is initially correlated with the environment, the state of 
the system in the future will depend on these correlations. This is the memory 
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of a system, and it is a key feature of non-Markovian dynamics [13j. We will 
show here that we can find the dynamics of the memory from the M-map. 

The memory of the system is obtained from the M-map by tracing over 
the system, indices r" and s". 

6r"s"y^r'r":s's" = Ure^r' aP%U*se,s' 0- (7-16) 



The last equation is the same as Eq. 3^ for the process map found by stochas- 
tic preparations. Thus, 

A. = 5^WMX',1,,,„. (7.17) 

This means, even though the M-map contains the correlated total state, we 
can find the process for the ideal stochastic preparation from it. 

Consider the following matrix composed of the matrices in Eqs. |7.10 
andlTlTl 

= c/ (5^(5,»,.mX^I,v' ) ($^'^-5rvMX'l,v' ) (7-19) 

V r" / \ rr' / 

= UreyaP% s" Paf3^se,s' /3 (7.20) 

e 

The last equation similar to the expression for the M-map, except the 
state of the system and the state of the environment are uncorrelated. Let us 
now write the total correlated state as 

P'' = P'®P' + X, (7.21) 
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where x is the correlation matrix [17J. Writing the M-map in terms of the 
correlation matrix we get 

^lv';s's" = Ure,r'a{p%s"Pal3 + Xr"a;s"/3)f^se,s'/3- (7.22) 

Now consider the following quantity: 

qrirs) _ T^irs) _ n (rs) /y r,q\ 

•^r'r";s's" ~ ■'^^r'r";s's" ^r'r";s's" yi-^O) 

Ure,r'a{Pr"s"Pal3 + Xr"a;s"l3 " Pr" s" P a(i)U*se,s' (5 (7.24) 

€ 

= ^ Ure,r'aXr"a;s"l3U*^,s,i3. (7-25) 
e 

The last equation gives us the reduced dynamics of the correlations between 
the system and the environment. If the system and the environment were 
initially uncorrelated, then % will be zero, i.e. no memory is present. We call 
% matrix the memory matrix. 

This is an important result for studying non-Markovian systems. It 
is important to keep in mind that, the memory matrix will be different for 
different preparations. But, if we want to study how the correlations act as 
memory in the uninterrupted system then we can just act with the identity 



preparation, defined in Eq. |7.11| as, 

X^{t) = %l.f^„.g,g/i'J'[.}^„.,,,^„ (7.26) 

= ^ UreyaXr'a,s'f3U*^^s'l3y (7.27) 

where x^it) = Tr£[f/x?7^]. The quantity above is the memory for the uninter- 
rupted system due to initial correlations with the environment. 
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7.3 An example with M-map and the memory matrix 



Now, let us consider an example of the M-map and the memory matrix. 



Let us use the same initial states (2.36) and the unitary operator (2.37) as 



before. The M-map is given by Eq. 7.4 and the memory matrix is given by 



Eq. 7.25 



If we make the identity preparation on the M-map then we get 



MoT® = pS(t) 
1 



(7.28) 



{I + cos {2Lut) ajaj — C23 cos {2Lut) sin {2Lut) o"i}, (7.29) 



which is the same as the final state of the system in Chap. [2]Eq. 2.40 While 
making the identity preparation on memory matrix we get 



1 



-{— C23 cos {2ujt) sin {2iut) ai}, 



(7.30) 
(7.31) 



where x is the correlation matrix defined in Eq. 7.21 and x i^) is the reduced 



part of the time evolved correlation matrix. This is precisely the contribution 
to due to the correlations with the environment. 
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7.4 Discussion 

The M-map and the memory matrix are difficult to obtain in practice. 
We have given a recipe to find M-map partially for a qubit in Apps. [B] and 
[C| Our recipe is based on projective preparations only. This issue needs to be 
further investigated using techniques such as direct characterization of quan- 
tum dynamics [3H1 [39]. Yet, the M-map and memory matrix % are powerful 
theoretical tools in their own right. The M-map allows us to determine var- 
ious output states for any preparation procedure, while the memory matrix 
% provides a quantitative definition for the memory effect due to the initial 
correlations for any preparation procedures. This allows us to investigate and 
optimize preparation procedures. 
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Chapter 8 



Analysis of experiments 



In this chapter we analyze two quantum process tomography experi- 
ments performed by Myrskog et al.[lB] and M. Howard et. al SH]. Our 
critiques emphasizes the importance of having a consistent procedure of state 
preparation. 

8.1 Experiment by Myrskog et al. 

In this experiment, quantum process tomography of the motional states 
of trapped ^^Rb atoms in the potential wells of a one dimensional optical lat- 
tice is performed. Only two bound bands are considered, which are labeled 
as states |0) and The states are prepared stochastically. An initial state 
of the system is the ground state \g) {g\, and from it states |r) (r|, \i) and 
the fully mixed state |l are prepared. The states |r) and \i) stand for the real 
and imaginary coherence states, which are prepared by applying appropriate 
unitary transformations on the ground state. This is achieved by displacing 
the lattice for the real coherence state and for the imaginary coherence state, a 
quarter-period delay is added after the displacement. The identity state is pre- 
pared by letting a superposition state decohere. The input states are allowed 
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to evolve and the output states are determine by quantum state tomography. 

The process map is found following the usual procedure laid out in 
Chap. |3j Since there are particles lost to the neighboring cells, the map is 
not required to be trace preserving. Based on this loss they also argue that 
the map can pick non-physical behavior (not-completely positive). The map is 
forced to be "physical" (completely positive) by using the maximum likelihood 
method |58j . 

In our terminology, the states prepared are P^^'~^\ P*^^'+\ P^^~\ and 
The three projective states are prepared by a single consistent stochastic 
preparation, while the fully mixed state is prepared by letting the state p(^'+) 
decohere. Which means an additional pin map is applied to prepare one of 



the input states. As we saw in the example in Sec. 5.2.2, this can lead to a 
not-completely positive and non-linear process map. 

Furthermore, the input states have varying values for polarization. 
Their data is listed in the table below. 





Pg 


Pi 


Pr 


Pi 


p(3,-) 


0.90 


0.60 


0.69 


0.69 


p(3,+) 


0.10 


0.40 


0.31 


0.31 


p(l>+) 


0.82 


0.59 


0.85 


0.63 


p(2,+) 


0.84 


0.58 


0.64 


0.37 



where pj are the experimentally prepared states projected onto projectors 
p(m) rjj^g polarization of the imaginary state in the a2 direction is very 
low. This could mean that it is correlated with the environment, while the 
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polarization of the ground state along the negative cxs direction is almost unity, 
meaning it is only weakly correlated at best. 

Here, we have two potential causes for the process map to have negative 
eigenvalues. The first problem is with the experimental procedure; applying 
multiple stochastic preparations that may affect the state of the environment 
differently. The second problem may be unavoidable; since it is extremely diffi- 
cult to prepare pure states in a setup like this. Even if the negative eigenvalues 
are due to the initial correlations, we do not have a prescription to obtain a 
process map in a consistent fashion. For the second example shown in Chap. 
[6} we assumed that the initial correlations were constant throughout the pro- 
cess. This does not seem to be the case here. Since the ground state is almost 
pure, while the imaginary state is clearly not pure. Therefore the correlation 
with environment for these two inputs must be different. This concludes our 
analysis of this experiment. Let us now consider the second experiment. 

8.2 Experiment by Howard et al. 

In this experiment, the system that is studied is an electron configura- 
tion formed in a nitrogen vacancy defect in a diamond lattice. The quantum 
state of the system is given by a spin triplet (5* = 1). Again we will write the 
initial state of the system and environment as p^^. 

The system is prepared by optical pumping, which results in a strong 
spin polarization. The state of the system is said to have a 70% chance of being 
in a pure state Or more mathematically, the probability of obtaining |0) 
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is Tr[|0) (01 p§^] = 0.7. 

Since the population probability is high, an assumption was made that 
the state of the system can be simply approximated as a pure state |0) 
From this initial state, different input states can be prepared by suitably ap- 
plying microwave pulses resonant with the transition levels. After prepara- 
tion, the system is allowed to evolve, and the output states are determined by 
quantum state tomography. With the knowledge of the input state and the 
measured output states, a linear process map is constructed. 

It was found that the hnear process map has negative eigenvalues, so 
the map was "corrected" using a least squares fit between the experimentally 
determined map and a theoretical map based on Hermitian parametrization 
[5U] . while enforcing complete positivity. 

However, if we do not regard the negative eigenvalues of the map as 
aberrations, then we should consider the assumptions about the preparation 
of the system more carefully. The assumption about the initial state of the 
system is: 

pS^^|0)(0|®r. (8.1) 

This is in effect a pin map. Along with the pin map, the stochastic transfor- 
mations are applied on the initial state to prepare the various input states; 



this is identical to the stochastic preparation method discussed in Sec. 4.3.1.1 

It is clear that the pure initial state assumption is unreasonable given 
our knowledge now of how the process is sensitive to the initial correlations 
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between the system and the environment. In effect the action of the pin map 
in this experiment is not perfect, and the pin map can be ignored. Then the 
process equation is: 



(8.2) 



where Q^"*) is the stochastic transformation that prepares the mth input state. 
In this experiment, Q^"^) is nothing more than a unitary transformation 
satisfying V^"^^ |0) = \^^"^^), where \ip^"^^) is the desired pure mth input state. 

We can write the unitary transformation for a two-level system as: 



(8.3) 



where (V'^^'^IV'J"^) = (0|0±) = 0. This defines as a transformation from 
the basis to the basis {|^|"^)}- The equation for the process becomes: 

+Tr£ \U IV^J"^) {ip^'^^l 



Therefore, since |0) = 0.7 to first approximation, the process is a linear 

mapping on the states IV'^"*-*) (-0^™^!- However, it is clear that if all terms 
are included, the process is not truly linear in the states jz/^^"*)) The 
negative eigenvalues are therefore a result of fitting results into a linear map 
when the process is not truly represented by a linear map. 
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Chapter 9 
Conclusion and future directions 



Our study of quantum process tomography started by noting that the 
dynamical map acting on a system can have negative eigenvalues. The dy- 
namical map has negative eigenvalues when the system is initially correlated 
with the environment. In the course of our studies of quantum process tomog- 
raphy, we showed that the preparation procedure cannot be neglected for any 
quantum system that interacts with an environment. These are the two major 
themes discussed in this dissertation. Though, along the way, we presented a 
method of quantum process tomography that is independent of the prepara- 
tion procedure. The map arising from this procedure lead us to an expression 
that quantifies the memory effect on the dynamics of the system due to the 
initial correlations with the environment. Determining the memory effect is 
an important task in coherence control. 

9.1 Future directions 

Though this dissertation has come to an end, the work goes unfinished. 
There are many open problems that arise from our work. We pointed out sev- 
eral causes for errors in quantum process tomography. We have not developed 
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a complete method of distinguishing and correcting these errors. This is an 
important problem that has not been fully addressed. 

Another problem that needs consideration is fully determining M-map. 
This is a challenging task, though there are some possibilities in resolving this 
problem by considering a scheme like the direct characterization of quantum 
dynamics. In this method a map is found with the aid of an ancillary system, 
which allows for direct estimation of the elements of a map without the aid of 
quantum state tomography. 

Another avenue of research that may prove fruitful is studying various 
preparation procedures using M-map. Furthermore, the memory matrix can 
aid in better better understanding of the various preparation procedures. This 
will especially be the case when the system is multi-partite. 
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Appendix A 



Quantum state tomography 



Quantum state tomography [T] is the procedure to experimentally de- 
termine the state of a quantum system. Suppose we want to determine the 
density matrix Q. To fully determine the state, we need to measure the mag- 
nitude of the polarization along the principal directions. For a qubit we need 
to make measurements along the cr^., cXy, and cr^ directions. 

However, a single experiment will not tell us the magnitude along any 
direction. Therefore we make many measurements of the magnitude along 
each of these direction. In general, the outcome of a set of the experiment 
along jth direction is written as 



which is the expectation value for the observable aj. Where cxj are the gener- 
alized Pauli-Gell-Mann-Tilma [63j matrices. Finally the total density matrix 
is given as follows 



Q 




(A.l) 
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Appendix B 



M-map process tomography 



We now develop a tomography procedure to determine the M-map. We 
will need to figure out a finite set of input states p(™) and the corresponding 
output states Q^™) that will allow us to determine the M-map. After-all, this is 
the objective of tomography; by performing measurements on a small number 
of select input states a complete description of a process can be obtained 
and predict the output state for any input state. The process equation for 
projective preparation^ takes the following form 



Q{m) _ _J_yr{rs) rp(m) /g , x 

^ ~ ^(m) r'r";s's"^r'r";s's" K'-'--^) 



~ ^(m) Vr";s's"-'s's" • \^-'^) 

The M-map is a large {d^ x d^) matrix, to make it more manageable, 



Eq. |7.3| can be interpreted as: 

(P(")| M |P(")) = r("^)g(™), (B.3) 
where P^™-) is now treated as a vector. In this form, M-map is a c?^ x d? matrix. 



where each element is a d x c/ matrix. We will call Eq. B.3 the ^'M-map process 
equation'. 



^Our procedure here only deals with projective preparations. 
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Since M-map is a Hermitian matrix, it has + (P) independent 



elements. Therefore ^(ci^ + (P) independent equations in the form of Eq. B.3 
are necessary to fully determine M-map. It is clear that neither an orthonormal 
set of d nor a linearly independent set of d"^ input states would provide sufficient 
equations to resolve the elements of the M-map. 

B.l The qubit case 

For a qubit system the projections, P^"^\ can be written in terms of 
three real parameter aj and Pauli spin matrices aj: 

P(-) = ^ (l + 4™ V,) . (B.4) 

Since P^*") is a projection then there is have the additional constrain ('^j'"'')^ " 
1. The matrices I and aj together forms a vector basis for this space. There- 



fore Eq. |B.4| is simply a vector decomposition of the projections in a fixed 

gives: 



basis. Taking this form for P^™) and substituting into Eq. B.3 



4^MqM ^ 4(P("*)|M|P('")) 

= (I|M|I) + aJ.'") {I\M\aj) + (afc|M|I) 

+4™)af)(a,|M|a,). (B.5) 

Observe that the terms (I|M|I), (I|M|o-j), (o-j|M|I) and {aj\M\ak), are simply 
the matrix elements of M-map in {I, aj} basis. We just need to find a set of 
projections P^*") that will allow us to solve for these matrix elements. 
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Consider the following specific projections defined as P^^'^^ = |(I±(jj) 
with j = {1,2, 3}. 

= (I|M|I) + {(Jj\M\aj) ± ((7,-|M|I) + {I\M\aj) . (B.6) 

Simultaneously solving the (+) and the (— ) equations above to gives the fol- 
lowing unknowns: 

(I|M|I) + {aj\M\aj) = 2 (r(^'+)Q(^'+) + r^^'-^Q^^'")) 
(I|M|(7,-) + ((7,-|M|I) = 2 (r^^'+)Q(^'+) - r(^'-)g(^'-)) . 

To obtain the cross terms {aj\M\ak) consider projections such as pO+'=+i'+) — 
+ -^aj + ^o^fc) for k> j which give: 

^{j+k+l,+)Q{j+k+l,+) ^ ^p(j+fc+i,+)| ^ |p(j+fc+i,+)^ 

= lim\I) + {aj\M\aj)) 
+^(+ (I|M|I) + K|M|afc)) 
+^((I|M|c7,) + (a,|M|I)) 

+^(+(Ma,) + (a,|M|I)) 
+^((a,|M|afe) + (a,|M|a,)). (B.7) 
Substitute the known terms and solve for the desired cross terms, 

{aj\M\ak) + {ak\M\aj) = -2 (l + V2j (r^^'+^Q(^'+^ - 2 (l + V2^ + r(^'+)Q('='+)) 

-2 (l - ^) (r(^'-)g(^'-) - 2 (l - x/2) + r('='-)g('='-)) 
_^8^(i+fe+i,+)g(i+fe+i.+) _ 
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In summary using the following nine projections, 

P(-+) = i (I + a,) , P(--) = ^ (I - a,) , P(^'+) = l[^ + ^-i + -^-2) 
P^''^^ = 1(1+ + -^a,) , P(«'+) = ^ f I + + 



and solving them simultaneously yields all desired matrix elements: (I|M|I) + 
{aj\M\aj), {I\M\aj) + {crj\M\I), and {aj\M\cxk) + {cxk\M\cXj) . 

This is not enough to fully determine the M-map, but these elements 
are sufficient to determine the output state for any input state. Using the 

can be rewritten as: 



property ^ -'- 



B.5 



j 

+ 5^4-)((I|M|a,) + (a,|M|I)) (B.8) 
j 

+ J]4-)4™)((a,|M|a,) + (a,|M|a,)). 

k>j 

Observe that the sums of the cross terms {aj\3Vl\ak) + (crfc|M|crj) can 
appear together because the coefficients a^J^^ are real. Also, the element (I|M|I) 
can always be paired with a diagonal element ((Tj|M|(Tj) as long as the state 
is a projection satisfying ^j{a^f^^)'^ = 1. The diagonal element (I|M|I) only 
has to be known if the system can be prepared directly to a mixed state such 
that E^r^f < 1- This may be accomplished by a generalized preparations 
If generalized preparations are allowed then just one more input state 
is needed, for example ^(I + ^cn), which gives another independent equation 
that can be solved to obtain (I|M|I). 
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Therefore the elements of the M-map found in Eq. |B.8 are all that 



are needed to describe the process. By measuring the outputs for the nine 
specified input states, the matrix M-map can be calculated. We now have a 
good quantum process tomography procedure for an open qubit system. 

B.2 Beyond one qubit 

Note that the choice of these nine states used above is not unique. 
The recipe which used to derive these nine states can be used in principle 
to derive other choices, and can also be partly generalized to d-level systems 
[63t However, there are some non-trivial obstacles to overcome for the 
generalization to d-\eve\ systems. In place of the Pauli matrices for two-level 
systems, the generalized Pauli-Gell-Mann-Tilma ^\ Hermitian traceless dxd 
matrices can be used |63l |65] to decompose the dxd density matrix, and 
this decomposition will also have only real coefficients. This trick eliminates 
certain degrees of freedom in the M-map that is otherwise difficult to deal 
with. Unfortunately, for d > 2, these real coefficients no longer satisfy just the 
simple constraint Ylji^^^j^^Y — 1- The additional constraints on the coefficients 
complicate the task of constructing the projections needed to simultaneously 
span the matrix elements of M-map. 
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Appendix C 



Experimental recipe for M-map 



Although we have estabhshed the M-map process tomography, let us 
make the ideas more concrete by developing a complete recipe for an exper- 
iment that can be used to determine whether a process is linear or given by 
the M-map. We will also show specifically how the M-map or linear map can 
be calculated from the measurement results. 

For the M-map process tomography nine input states are necessary. 
For the nine states derived in App. [B] the first six states are three pairs of 
orthonormal projections, but the last three are not. If we use projections 
to prepare the states, we will need use them as given by orthonormal pairs. 



Therefore, let us instead use twelve projections, nine from Eq. B.8 and three 
orthogonal to the last three in that equation: 

p(6,-) = 

These twelve projections are neatly grouped into six different sets of 
orthonormal pairs. If the states are prepared using projections, these would 

















7=2-') 
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correspond to measurements in the ±cri, ±(T2, ±cr3, cri ± (T2, cxi ± (J3 and 
(J2 ± (T3 directions. Having twelve states is more than is necessary for the JVC- 
map process tomography, but the extra states and the corresponding output 

states can serve as consistency checks. 

After recording the corresponding output states for all twelve input 
states, we can verify if the process is linear by checking the following eight 
linear sum rules (that the input states satisfy) have to be satisfied by the 
corresponding output states: 

Q{3,+) + Q3,-) ^ Q{k,+) ^ Q(k,-) (fQj, k^l,2, 3) 

If the eight sum rules are satisfied, then the process is not described by 
the M-map, and we can be confident the process is described by a linear map. 
The linear map A can then be computed as follows: 

a(^(i + p.<^.)) = {Q^''^^ + Q^''-^) + Pj {Q^'^^^ - Q^''-^) (C.2) 

If the eight sum rules are not satisfied, then the process may be given 
by the M-map. We will attempt to verify that the process given by the M-map 
and calculate it. 
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Note that the probabihties r^"^^ = Tr[p^^p(™')] associated with each 
preparation should be found experimentally. The probabilities should be com- 
plete for an orthonormal set of projections, in other words: r^^'^'> + r^^~^ = 1. 
Therefore the probabilities can be calculated from the fraction of the (+) states 
as compared to the (— ) states, for all preparations made in the same direction. 

To be certain that the process is given by the M-map, we can check 
if the three additional states we included evolve in a way that is consistent 
with the M-map derived from the other nine states. The following equations 
are derived from this condition. If these equations are satisfied, we can be 
confident that the process is given by the M-map: 

V2 (r(^'+)g(^'+) - r(5'-)g(5'-)) = r(i'+)g(i'+) - r(i'-)g(i'-) 

+^(3,+)g(3,+) _ ^{3,-)g(3,-) 
V2 (r(^'+)g(^'+) - r^'^'-^g^^'-)) = r(2.+)g(2.+) _ r(2'-)g(2.-) 

+^(3,+)g(3,+)_^(3,-)g(3,-)^ 

If the conditions above are satisfied then the process is given by the M-map, 
then the M-map can be computed by following the recipe in App. [B} 

Once the matrix elements of M-map are determined, the evolution of 
a generic state X = |(I + J2j ^j^j) is given by: 

4(X|M|X) = r(^)gW (C.3) 
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= ^x|((I|M|I) + (a,|M|a,)) (C.4) 
j 

+ 5^x,((I|M|(T,) + ((Tj|M|I)) (C.5) 
j 

+ J2^j^k{{(rj\M\ak) + {ak\M\a,)). (C.6) 
k>j 

Note that since Q is a normalized state, the normahzation constant r 
is the measurement probabihty r = Tr[Xp^^]. Ahhough we had not exphcitly 
mentioned this before, the M-map contains all information about the mea- 
surement probabilities, that is why we needed the measurement probabilities 
r^"^'^ to calculate the M-map. 

Finally, note that if both the test for linearity and the M-map fail, 
then the process cannot be consistently described by either a linear map or 
the M-map. The experiment then should be carefully analyzed for problems 
such as any non-linear dependence that may have been introduced if the input 
states are not accurately prepared or if there is some dependence of p^^ on the 
prepared state. 

If the process is given by a linear process map, then the nine input states 



(Eq. B.8) are over complete; only four input states are needed to determine 
a linear process map. This discrepancy is summarized by the following linear 
sum rules: 

p(l.+) + p{l -) = p(2,+) _^ p(2-) ^ p(3,+) _^ p(3-)^ 
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p(5,+) ^ 




p(6,+) ^ (1 L ) /'p(l,+) + p(l-)) + J_ /'p(2,+) I p(4,+)) 

V2 V2y^ ^ 

If the process is linear, then the output states must satisfy the same 
sum rules, which are obtained from the above equations by suitably writing Q 
in place of P. If these sum rules are not satisfied, then the process is not linear. 
However satisfying the sum rules is necessary but not sufficient to determine 
if the process is hnear; the M-map can still be constructed from this set of 
input and output states without any contradictions. Therefore an additional 
input state, distinct from the above nine input states should be tested. 
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Appendix D 



Notation 



a : 


: Bloch vector 


A : 


: stochastic map 


!B : 


: dynamical map 


Cjk 


: correlation parameters 


X - 


correlation matrix 


H : 


: Hamiltonian 


I : 


: indentity matrix 


3 : 


: identity map 


X : 


: memory matrix 


A : 


: linear process map 


M : 


: Kuah's dynamical M-map 


n^"") : 


: mth local rotaion 


rp{m) , 


: mth preparation map 


p{m) , 


: mth input states 




: mth output states 
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p : state of the ancillary system 

: state of the environment 

p^ : state of the system 

p^^ : state of the system plus environment 

Uj : Pauli spin matrices 

© : pin map 

U : unitary operators 

X : an arbitrary input state 
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